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We review various field theory approaches to the description of neutrino oscillations in vacuum 
and external fields. First we discuss a relativistic quantum mechanics based approach which involves 
the temporal evolution of massive neutrinos. To describe the dynamics of the neutrinos system we 
use exact solutions of wave equations in presence of an external field. It allows one to exactly 
take into account both the characteristics of neutrinos and the properties of an external field. In 
particular, we examine flavor oscillations an vacuum and in background matter as well as spin flavor 
oscillations in matter under the influence of an external electromagnetic field. Moreover we consider 
the situation of hypothetical nonstandard neutrino interactions with background fermions. In the 
case of ultrarelativistic particles we reproduce an effective Hamiltonian which is used in the standard 
quantum mechanical approach for the description of neutrino oscillations. The corrections to the 
quantum mechanical Hamiltonian are also discussed. Note that within the relativistic quantum 
mechanics method one can study the evolution of both Dirac and Majorana neutrinos. We also 
consider several applications of this formalism to the description of oscillations of astrophysical 
neutrinos emitted by a supernova and compare the behavior of Dirac and Majorana neutrinos. Then 
we study a spatial evolution of mixed massive neutrinos emitted by classical sources. This method 
seems to be more realistic since it predicts neutrino oscillations in space. Besides oscillations among 
different neutrino flavors, we also study transitions between particle and antiparticle states. Finally 
we use the quantum field theory method, which involves virtual neutrinos propagating between 
production and detection points, to describe particle-antiparticle transitions of Majorana neutrinos 
in presence of background matter. 
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Neutrino physics is one of the most rapidly developing 
area of high energy physics especially after the great sue- 
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cess in experimental studies of neutrino properties. In 
the first place we should mention the investigation of as- 
trophysical neutrinos since they play an important role in 
the evolution of various astronomical objects like stars, 
supernovas, quasars etc. In the course of recent experi- 
ments for the detection of solar neutrinos [l| it was re- 
vealed that transitions between different neutrino flavors, 
neutrino oscillations, are the most plausible theoretical 
explanation of the solar neutrinos deficit Q ■ Flavor oscil- 
lations were also observed in the experimental studies of 
atmospheric neutrinos [H . It is worth noticing that there 
are numerous attempts to detect neutrinos from outside 
the solar system (see, e.g., Ref. [4|), however, presently 
only SN1987A supernova neutrinos were observed @. 

Besides natural sources, one can use neutrinos pro- 
duced by an accelerator or a nuclear reactor to study 
oscillations of these particles Q . In this case we control 
the flavor content of both initial and final fluxes, i.e. it 
is the best strategy to examine neutrino oscillations. Be- 
sides the studies of neutrino oscillations, some accelerator 
based experiments are dedicated to the investigation of 
neutrino interactions (see, e.g., Ref. Q). 

The existence of neutrino oscillations is a direct indica- 
tion to the facts that neutrinos are massive particles and 
there is a mixing between different neutrino generations. 
There are multiple theoretical mechanisms to generate 
masses and mixing of neutrinos Q to fit the data of the 
aforementioned neutrino oscillations experiments. 

It was found that besides the possibility of flavor con- 
versions in vacuum, various external fields, like inter- 
action with background matter [9] or with an external 
magnetic field [10], can significantly influence neutrino 
oscillations. For example, the resonant enhancement 
of neutrino oscillations in matter, Mikheev-Smirnov- 
Wolfcnstein (MSW) effect, plays an important role in 
the solution of the solar neutrino problem [ll| . It should 
be noted that the standard model of electroweak in- 
teractions does not imply the mixing between different 
neutrino flavors when particles propagate in background 
matter. However the possibility of nonstandard neutrino 
interactions, which can cause transitions between neu- 
trino flavors even at the absence of vacuum mixing, was 
recently discussed [12 1. 

Since neutrinos are unlikely to have a nonzero elec- 
tric charge [l3| , the interaction with an external electro- 
magnetic field may be implemented due to the presence 
of anomalous magnetic moments. Note that, unlike an 
electron which has a vacuum magnetic moment, electro- 
magnetic moments of a neutrino always arise from the 
loop corrections. The contributions to neutrino electro- 
magnetic moments in various extensions of the standard 
model are reviewed in Ref. flij ]. 

It should be noted that neutrino electromagnetic prop- 
erties have rather complicated structure. First, in the 
system of mixed neutrinos there can be both diagonal 
and transition magnetic moments. In presence of an ex- 
ternal electromagnetic field the former are responsible 
for the helicity flip within one neutrino generation and 



the latter cause the change of neutrino flavor (spin fla- 
vor oscillations). Second, electromagnetic properties of 
Dirac and Majorana neutrinos are completely different. 
In general case Dirac neutrinos can have all kinds of mag- 
netic moments, whereas Majorana neutrinos possess only 
transition magnetic moments which are antisymmetric 
in neutrino flavors 15]. Nowadays there is no univer- 
sally recognized experimental confirmation of the nature 
of neutrinos [l6| . Thus it is important to study the evolu- 
tion of both Dirac and Majorana neutrinos in an external 
electromagnetic field. 

The indication to the Majorana nature of neutrinos 
should be experimental confirmation of the existence of 
neutrinoless double /3-decay (0^2/3), which is the mani- 
festation of oscillations between neutrinos and antineutri- 
nos [TtJ ■ This kind of transitions is possible only if neutri- 
nos are Majorana particles. Despite numerous attempts 
to detect (0^2/3)-decay in a laboratory were made [l8| . 
no confirmed events are known presently. 

The great variety of evidences of neutrino oscillations 
requires a rigorous theoretical explanation of this phe- 
nomenon. The conventional approach to the description 
of neutrino oscillations is based on the quantum mechan- 
ical evolution of neutrino flavor eigenstates governed by 
an effective Hamiltonian [l9j]. This intuitive approach 
may be easily extended to the description of neutrino 
oscillations in background matter Q and spin flavor os- 
cillations in an external magnetic field [20] giving one 
a reasonable description of neutrino oscillations in these 
external fields. 

Neutrinos are supposed to propagate as plane waves 
in the quantum mechanical approach. However, due to 
uncertainties in production and detection processes, neu- 
trinos seem to have some distribution of their momenta, 
i.e. they propagate like wave packets [2l|. The Gaus- 
sian distribution of the neutrino momentum is typically 
discussed (22[. However the actual form of the distribu- 
tion strongly depends on the production and detection 
processes (see, e.g., Ref. 23]). The necessity of the wave 
packets treatment of neutrino oscillations is discussed in 
Ref. [24|. The details of this approach including its ex- 
tension to the Dirac theory were considered in Ref. [25| . 

The attempts to reproduce the quantum mechani- 
cal transition probability formula in vacuum were made 
in Refs. [26], [27| treating massive neutrinos as virtual 
particles propagating between production and detection 
points. Using this quantum field theory method one deals 
with the observables of charged leptons rather than with 
the characteristics of mixed neutrinos. The special Gaus- 
sian form of the source and detector within this quantum 
field theory approach was discussed in Ref. 28]. The 
analysis of relativistic effects in neutrino oscillations us- 
ing this approach was made in Ref. [29| . 

In Ref. [30[ neutrino flavor oscillations in vacuum are 
described on the basis of time evolution of the Fock states 
of flavor neutrinos. The time dependent transition prob- 
ability formula, obtained within the S'-matrix approach, 
was studied in Ref. [311 ] . The quantum mechanics analy- 
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sis of the neutrino detection process and the possibility to 
obtain the nonstandard transition probability formulas is 
considered in Ref. [32]. 

We contributed to the field theoretical description of 
neutrino oscillations in Refs. [33ll40l |. We developed an 
approach based on relativistic quantum mechanics and 
applied it for the description of neutrino flavor oscilla- 
tions in vacuum [33[ , in background matter [34|, [36| , and 
to spin flavor oscillations in an external electromagnetic 
field [H, Hi], [H 53]. This method is based on the first 
quantized neutrino wave packets which was also studied 
in Ref. (4l| to study neutrino oscillations in vacuum and 
in an external electromagnetic field. 

In frames of the relativistic quantum mechanics 
method we formulate the initial condition problem for 
the system of flavor neutrinos and study the subsequent 
time evolution of neutrino wave packets. When we dis- 
cuss neutrino propagation in an external field, we use 
an exact solution of the corresponding wave equation in 
presence of this external field. With help of this method 
we could reproduce the Schrodinger type evolution equa- 
tion, which is used in the standard quantum mechanical 
description of Dirac and Majorana neutrino oscillations, 
and discuss the correction to the standard approach. 

In Ref. [38| we studied neutrino flavor oscillation in 
vacuum within the external classical sources method. 
Note that the analogous approach to the description 
of neutrino oscillations was also discussed in Ref. [42j . 
Within this approach we studied the spatial evolution 
of flavor neutrino waves emitted by external classical 
sources. We could describe the evolution of both Dirac 
and Majorana neutrinos in vacuum. 

In the present work we review our recent achieve- 
ments in theoretical description of neutrino flavor oscil- 
lations. This work is organized as follows. In Sees. HlWl 
we discuss oscillations of Dirac neutrinos in frames of 
relativistic quantum mechanics method. We start with 
the description of neutrino flavor oscillations in vacuum 
(Sec. |TT|. Then we consider oscillations in background 
matter (Sec. IIIip . where we also study oscillations in 
case of nonstandard neutrino interactions with matter. 
In Sec. IIVI we apply the relativistic quantum mechanics 
method to the description of spin flavor oscillations of 
Dirac neutrinos in an external magnetic field. Finally, 
in Sec. |V] we discuss the most general situation of spin 
flavor oscillations in matter under the influence of an ex- 
ternal magnetic field. In Secs. lVII and lVIII we examine ap- 
plications of the relativistic quantum mechanics method 
to the description of propagation and oscillations of as- 
trophysical neutrinos in a magnetized envelope after the 
supernova explosion. 

In Sees. I Villi and IPXl we analyze the evolution of mas- 
sive mixed Majorana neutrinos in vacuum as well as 
in matter and magnetic field in frames of the relativis- 
tic quantum mechanics approach. Besides oscillations 
among different neutrino flavors we also consider transi- 
tions between neutrino and antineutrino states which are 
allowed if neutrinos are Majorana particles. In Sec. [X] 



we apply the general results to the studies of oscilla- 
tions of astrophysical neutrinos in the supernova enve- 
lope supposing that neutrinos are Majorana particles and 
compare them with the Dirac neutrino case studied in 
Sec. EH 

Then, in Sec. IXI[ we formulate an alternative for- 
malism for the description of neutrino flavor oscillations 
which is based on the mixed massive neutrinos emission 
by classical sources. We examine the spatial distribu- 
tion of flavor neutrino wave functions in case of localized 
sources. Note that vacuum oscillations of both Dirac 
fSec. lXI Al and Majorana fSec. lXIB]) neutrinos are stud- 
ied. When we discuss Majorana neutrinos case, we also 
consider the possibility of transitions between particles 
and ant ip articles. 

Finally, in Sec. I XIII we consider the transitions among 
neutrinos and antineutrinos in frames of the quantum 
field theory, treating neutrinos as virtual particles prop- 
agating between macroscopically separated production 
and detection points (see also Refs. (26l - [29| ). In par- 
ticular we are interested in the influence of background 
matter on the oscillations process since this kind of 
transitions typically happens inside a nucleus (see, e.g., 
Ref. [17]) and one cannot neglect the presence of dense 
nuclear matter in neutrino oscillations. 

In Sec. IXIIII we summarize our results. Several tech- 
nical issues are considered in Appendices IaTIdI in order 
not to encumber the description of the main results in 

sees, imixm 



II. DIRAC NEUTRINOS IN VACUUM 

In this section we discuss the evolution of mixed flavor 
neutrinos in vacuum, i.e. at the absence of external fields, 
using the relativistic quantum mechanics approach [33j j . 
First we formulate the initial condition problem for the 
system of flavor neutrinos. We exactly solve this prob- 
lem and find the time dependent wave functions and the 
transition probability. We also discuss the validity of the 
developed formalism. 

Without loss of generality we study the situation of 
only two particles v = (v a , vp), where a and (3 can stay 
for electron, muon or T-neutrinos. In various theoretical 
models (see, e.g., Ref. [43j]) bigger number of flavor neu- 
trino fields N v > 3 is proposed. The case of arbitrary 
number of neutrino fields can be considered in our for- 
malism by straightforward generalization. The Lorentz 
invariant Lagrangian for this system has the following 
form: 

C= ^ ^xvy^dnVx- ^2 m AA' v\v\i , (2.1) 

\=a,P \X'=a,/3 

where 7^ = (7 ,7) are Dirac matrices and (m\\>) is the 
mass matrix which is not diagonal in the flavor neutrinos 
basis. The nondiagonal elements of this matrix m a p = 
mp a correspond to the mixing between different neutrino 
flavors. 
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To describe the time evolution of the system (|2.1[) we 
formulate the initial condition problem for flavor neutri- 
nos V\, 



^ A (r,i = 0) = ^ O) (r), \ = a,P 



(2.2) 



where (r) are known functions. Eq. (|2.2p means that 
initial field distributions of flavor neutrinos are known 
and we will search for their wave functions at subsequent 
moments of time: v\(v, t > 0). The situation when only 
one of the flavor neutrinos, e.g. belonging to the type 
"/3" , is present initially corresponds to a typical neutrino 
oscillations experiment: one looks for the initially absent 
neutrino flavor "a" in a beam consisting of neutrinos of 
the flavor "/3" . 

To solve the initial condition problem (|2.2I) for the sys- 
tem (I2.ip we introduce the mass eigenstates neutrinos ip a , 
a = 1,2, 



^2 ^AaV'a, 



(2.3) 



where the matrix (U\ a ) is chosen in such a way to di- 
agonalize the mass matrix (m\\>). The eigenvalues of 
the matrix (raw), which are real and positive, have the 
meaning of the masses of the fields ip a . We define them 
as m a . 

The Lagrangian formulated in terms of the flavor neu- 
trino fields does not provide any information about the 
nature of neutrinos, i.e. whether neutrinos are Dirac 
or Majorana particles, since in general case it is writ- 
ten using the two component left and right handed 
spinors (4J, [45j . Only when we introduce the mass eigen- 
states and analyze the structure of resulting mass matri- 
ces, we can reveal the nature of neutrinos. We suppose 
that in our case the fields ip a are Dirac particles. 

For the case of only two neutrinos system the mixing 
matrix (U\ a ) in Eq. (|2.3p has the form, 



(Uxa) 



cos u — sin ( 



(sin 9 



(2.4) 



where 9 is the vacuum mixing angle. 

The Lagrangian (|2.1I) expressed in terms of the fields 
ip a reads 



C = ipaij-^dp - m a )ip a . 



(2.5) 



o=l,2 



The Lagrangian (|2.5I) should be supplied with the initial 
condition, 

^ a M = 0)=^ o >(r), ^ 0) = (O^ 0) , (2-6) 

which follows from Eqs. ([2^21 and (f!Q|) . 
The Dirac equations, 



iip a = (ap + Pm a )ip a , 



(2.7) 



which result from the Lagrangian (|2.5I) . reveal that the 
fields ip a decouple. In Eq. (I2.7[) we use the standard def- 
initions for Dirac matrices a. = 7°7 and (3 = 7 . The 
solution of Eq. (|2.7[) can be found in the following way: 



d3 P P ipr 



(2tt) 3 / 2 
C=±i 



, (2.8) 



where E a — -\/|p| 2 + m 2 is the energy of a massive neu- 
trino in vacuum, £ = ±1 is the helicity of massive neutri- 
nos, and Ua and Va are the basis spinors corresponding 
to a definite helicity. 

In the relativistic quantum mechanics approach to the 
description of neutrinos evolution 33] the coefficients 

and ba are c-number quantities rather than operators 
acting in the Fock space. Our task is to find these coeffi- 
cients. Since the fields ip a are independent, the values of 
the coefficients ai and 60 depends only on the initial 
condition (|2.6[) . 

Using the orthonormality of the basis spinors, 



COUO =( y f#fi u (Ot v (0 =0j (2.9) 



= <>(,'(,: 

we find the coefficients a! and bffl in the form 



a (2tt) 3 / 2 a a 

&(0 = - v^ib^(v) 

(27r) 3 / 2 



^ 0) (P)= / cl^e-^^W, 



(2.10) 



(2.11) 



where 



is the Fourier transform of the initial condition (|2.6I) . 

With help of Eqs. (|2.8p - (|2.1ip we arrive to the expres- 
sion for the wave function of the neutrino mass eigen- 
states. 

Mr,t) = J d 3 r'S„(r' -r,t)(-i 7 >i 0) (r'), (2.12) 
where 

S„(r, t) = (yfdp + m a )D a (r, t), (2.13) 
is the Pauli-Jourdan function for a spinor particle and 



d 3 p ■ sin E a t 



E a 



(2.14) 



is the Pauli-Jourdan function for a scalar particle. 
It is convenient to rewrite Eq. (|2.12l) in the form, 



d 3 p 

(2tt) 3 



e i P r 5 a (-p,t)(-i 7 °)^ 0) (p), (2.15) 
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where 



C=±i 



(0 

cos E a t — i 



(C)tpiBa* 

sin _E a £ 

E a 



(ap + /3m Q ) 



(2.16) 



is the Fourier transform of the Pauli-Jourdan func- 
tion (|2.13f) . To derive Eq. (|2.16[) we use the summation 
over the helicity index formulas (46[, 



C=±i 



1 

2Ea 



(ap + f3m a ) 



(2.17) 



which are consistent with the normalization of the basis 
spinors (|2.9[) . 

Now let us specify the initial condition (|2.2[) . We sug- 
gest that initially very broad wave packet is present, i.e. 
the coordinate dependence of the wave functions u^\r) 
is close to a plane wave corresponding to the initial mo- 
mentum k. Moreover we choose the situation when only 
one flavor neutrino is present, 



^°)(r)=0, ^(r) = e ik ^°>(k), 
,(o) 



(2.18) 



where v a (k) is the coordinate independent normaliza- 



tion spinor, |iA°^(k)| 2 = 1 



Using Eqs. (|2"3)l . (|2~i|) . and (l2TTi|) - (l2TTg)) we obtain the 
wave function of the initially absent flavor neutrino v a as 



v a (r,t) =e lkr sin6'cos#^ cos E\t — cos E 2 t 



.smEit. . 
i — - — (ak + /3mi) 
Ei 

. sin E 2 t , 
l — — — 

E 2 



(ak + /?m 2 H^ 0) (k), (2.19) 



where the energies are the functions of the initial mo- 
mentum, E a = ^/|k| 2 + m 2 . 

With help of Eq. (|2.19|) we get the transition probabil- 



ity P, 



P, 



(t) 
(*) 



\is a (r,t)\ 2 as 



= sin 2 (29) 
1 



sin 2 (<K) - sin($i) cos(<rt) 



2 
O 



where 



Ex 



2 

E\ + E-2 



k 2 

4 

k 4 



8m 



f sin E x t 



k 2 



smEjt 



(2.20) 



4fc 



(2.21) 



and Sm 2 



■ m\ . The quantity $ has the meaning of 



the phase of neutrino oscillations in vacuum. Note that 
the coordinate dependence is washed out from Eq. (|2.20j) 
since we study a very broad initial wave packet. 

In the untrarelativistic limit (I2.2ip the main term in 
Eq. (|2.20|) resembles the usual transition probability for 
neutrino oscillations in vacuum. The correction to the 
main result is suppressed by the factor m 2 /k 2 <C 1. Using 
Eq. (|2.21[) we can represent Eq. (|2.20[) in the following 
form: 



Pup^vS*) = sin 2 (20) 



Sm 



8m 2 
Ifc 2 



4A- 



t sin(2fct) 



(2.22) 



where we drop small terms ~ m\/k . 

Now the leading term in Eq. (|2.22j) reproduces the well 
known transition probability for neutrino oscillations in 
vacuum derived in frames of the quantum mechanical ap- 
proach [47[ . The correction to the leading term, which is 
a rapidly oscillating function on the frequency ~ fc, was 
first studied in Ref. [3(| in frames of the quantum field 
theory approach to neutrino oscillations. We obtained 
analogous result using the relativistic quantum mechan- 
ics method (see also Ref. [33[). This correction to the 
leading term in the transition probability results from 
the accurate account of the Lorentz invariance. 

Note that Eq. (|2.22|) is invariant under the mi O m 2 
transformation. It means that one cannot obtain the in- 
formation about neutrino mass hierarchy studying neu- 
trino oscillations in vacuum even taking into account the 
correction to the leading term in the transition probabil- 
ity. 

Now let us discuss the possibility of initial conditions 
which differ from the plane wave (|2 . 1 8|) . If the initial 
wave function is localized in a spatial region with a typ- 
ical size Lq and we measure a signal in the wave zone 
|r| Lq, the dependence on the particle masses m a 
is washed out from the Pauli-Jourdan functions (|2.14[) 
and (|2.13[) (see Refs. [H, IHjj). Thus neutrinos with spa- 
tially localized initial wave packets evolve in the wave 
zone like massless particles, which are known not to re- 
veal flavor oscillations. Therefore the initial wave packet 
should be sufficiently broad. 



III. DIRAC NEUTRINOS IN BACKGROUND 
MATTER 

In this section we use the formalism developed in 
Sec. [Til to study the evolution of the system of mixed fla- 
vor neutrinos v = (y a , vp) in background matter [34H36|. 
We formulate the initial condition problem for this sys- 
tem and solve it for ultrarelativistic neutrinos. The case 
of the standard model neutrino interactions is studied in 
details. Then we also analyze the dynamics of neutrino 
oscillations in presence of the nonstandard interactions 
which mix neutrino flavors. 
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The neutrino interaction with matter can be rep- 
resented in the form of an external axial-vector field 
f\y O IH3|- A s in Sec. HJ we start from the Lorentz 
invariant Lagrangian, 



\=a,/3 



Y ^(m X y + f^)^, (3.1) 

\,\'=a,l3 



i7°7 1 7 2 7 3 , and keep the 



where 7^ = 7/^(1 ~ 7 5 )A 7 5 
same notation for the mass matrix as in Sec. [TTJ 

Note that in general case the axial vector field f\y can 
be nondiagonal in the neutrino flavor basis. The appear- 
ance of nondiagonal elements (/„«) of this matrix is the 
indication to the presence of nonstandard neutrino inter- 
actions which mix neutrino flavors since in the standard 
model of electroweak interactions only diagonal elements 
of this matrix can appear. The time component of di- 
agonal elements of this matrix, is proportional to 
the density of background matter and spatial component, 
fx a, to the mean velocity and polarization of background 
fermions. The details of the averaging over the back- 
ground fermions are presented in Ref. [50| 

To describe the evolution of the system (|3.1[) we for- 
mulate the initial condition problem for flavor neutrinos, 
with the initial wave functions having an analogous form 
as in Eq. (|2.2[) . To solve the initial condition problem 
we introduce the set of neutrino mass eigenstates ifi a 
[see Eqs. (|2.3j) and (|2.4[) ] to diagonalize the mass matrix 
{m\x')- As in Sec. [II]we suggest that mass eigenstates 
tp a are Dirac particles. 

The Lagrangian (|3.1[) expressed in terms of the fields 
ip a has the form, 



0=1,2 



m a )if) a 



E 

a, 6=1, 2 



Sa5^7>, (3-2) 



where 



a&) = ^(/&')tf=($ (3-3) 

is the external axial-vector field expressed in the mass 
eigenstates basis. 

One can derive Dirac equations for the neutrino mass 
eigenstates, 

itpa =H a if>a + a, 6 = 1,2, a^t, 

Ua=ap + Pm a + Pfrg£, V = pftg*, (3.4) 

directly from the Lagrangian (|3 . 2[) . It should be noted 
that Dirac equations for different neutrino mass eigen- 
states are coupled because of the presence of the inter- 
action V. Studying an exact solution of the Dirac equa- 
tion for a massive neutrino in background matter we can 
exactly take into account the contribution of the term 
Plu9a m t° the dynamics of a particle. On the contrary, 
the term proportional to V, which mixes different mass 



eigenstates, should be studied perturbatively. Neverthe- 
less one can account for all terms in the perturbativc 
expansion for ultrarelativistic neutrinos. 

We will study the case of nonmoving and unpolarized 
matter which corresponds to g ao = 0. The matrix (g^ b ) 
has only time component now, 



(9ab) = (g^b) 



9i 9 

g 52 



(3.5) 



where we introduce the new notations, 



g = 9i2 



_ = 9aa and 

(721- If the background matter is nonmoving 
and unpolarized, the Hamiltonian ~H a commutes with the 
helicity operator (Sp)/|p|, where £ = 7°7 5 7, and we 
can classify the states of massive neutrinos with help of 
the eigenvalues of the helicity operator ( = ±1. 

The general solution of Eq. (13.41) can be presented in 
the following way: 



ip a (r,t) =e 



= p -\g a t/2 



d 3 p 



ipr 



(2tt) 3 / 2 

x Y, [aPW4 C) exp(-iS(°*) 
C=±i 

+ b®(t)vP exp(+LB(0t)], 



(3.6) 



where a„ and b a ^ are the undetermined nonoperator 



coefficients [see Eq. (|2.8p ]. which are, however, time de- 
pendent now because of the presence of the term V in 
Eq. (031). 

The energy spectrum Ea in Eq. (13.61) was found in 
Ref. 



4° = v(ipi-cw2r 



(3.7) 



for the case of nonmoving and unpolarized neutrinos. 
The basis spinors Ua and Va in Eq. (|3.6[) are the eigen- 
vectors of the helicity operator (Sp)/|p|, with the eigen- 
values £. As an example, we present here the basis 
spinors which correspond to an ultrarelativistic particle 
propagating along the z-axis, 





(3.8) 



where we omit the subscript "a" since we neglect the 
neutrino mass in Eq. (I3.8[) . Basis spinors corresponding 
to arbitrary energies were also found in the explicit form 
in Ref. [H|. 

Now we should specify the initial condition. We can 
choose it as in Eq. (I2.18[) . with k = (0, 0, k) and k ^> m a . 



It is also convenient to take v 



(0) 



[see Eq. 43U)] 
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Such an initial wave function corresponds to a neutrino 
propagating along the z-axis, with the spin directed op- 
posite to the particle momentum, i.e. it describes a left 
polarized neutrino. 

If we put the ansatz (|3.6|) in the wave equations (|3.4I) , 
we get the following ordinary differential equations for 
the functions a~a\t) and b^\t): 



v(0 = 



e i(ff.-m)t/a exp 



E K 



(iEPt) 



exp ( -iE^'h 



C'=±i 



+ bpv^ exp (lE^t) 



ibCO 



:e i(Sa-Sb)t/2 exp 



iEPt) v«»V 



E K 



C'=±i 



+ 6 



exp 



exp (i^ C ' ) t 



(3.9) 



To obtain Eq. (|3.9j) we use the orthonormality of the 
basis spinors p.8[) [see Eq. (|2.9p ]. We should supply the 
Eq. p. 91) with the initial condition, 



4°(o) 



(2tt)3/ 2 ^ ' 



cos 6* 
(27t) 3 / 2 



(3.10) 



that result from Eqs. (|2.6[) and (13.61) . If we study an 
arbitrary wave packet initial condition rather than the 
plane wave distribution for (r) (|2.18l) , we have to re- 
place in Eq. (|3.10[) with the Fourier transform of the 
initial wave function u^\r). 

Taking into account the fact that («'C|V|iA ') = 0, 
we get that equations for a^\t) and b£\t) decouple, 
i.e. the interaction V does not mix positive and negative 
energy eigenstates. In the following we will consider the 
evolution of only ai (i) since the dynamics of 6 a (i) is 
studied analogously. 

The only nonzero matrix elements of the potential V in 
Eq. (pT9|) are (u - |V|u"~) = (v+\V\v+) = g, which result 
from Eq. (|3.8I) . Finally Eq. (|3.9[) are reduced to the ordi- 
nary differential equations only for the functions a~ (t), 



--a b gexp[i{E a - E b + (g a - g b )/2}t], 
a, b = 1, 2, a =/= b. 



(3.11) 



It follows from Eq. p. Ill) that equations for the func- 
tions a~ and a+ (not shown here), decouple since the 
interaction with background matter conserves the parti- 
cle helicity and we have chosen the initial condition cor- 
responding to a left polarized particle. Indeed one can 
obtain from Eq. (|3.10|) that the functions a+ (i) are equal 
to zero at t = 0. 



The solution of Eq. (|3.11[) can be expressed in the form 
(see Appendix [A"]) , 



<z7/(t) =Fa7/(0) + 6^(0), 
a 2 -(t) =F*a 2 -(0)-G*ar(0), 



F = 
G = 



cos f2_£ — i- 



2f7_ 



■ sin fit 



exp (kj_i/2), 



sin£7_£exp {iu}—t/2), 



(3.12) 



where = \J g 2 + (w_/2) 2 and w_ = — i5 2 + (51 — 
92)/2. 

Using the identity (w+ <g> w +t ) i^ 0) = [sec Eq. (|3T8]) ] 
as well as Eqs. (|23]l . (|2~4| . p~6| . and (|3~T2|) we arrive to 
the wave function of the flavor neutrino v a , 



v a (z, t) — — iexp (— ier_t + ifcx 



sin Q _t 



x [p cos 26* + (cj_ /2) sin 20] i/£ 

■ o"" 



(o) 



(3.13) 



where <t_ = (£f + E^)/2 + (g x + g 2 )/2. Note that 
Eq. p,13p is the most general one which takes into ac- 
count the nonstandard interactions of relativistic neutri- 
nos with nonmoving and unpolarized matter of arbitrary 
density. 

Now let us discuss the standard model neutrino inter- 
actions with background matter. In this case the matrix 
{f XX') 1S diagonal: f\\i — /a^aa'- Since we study the 
nonmoving and unpolarized matter the spatial compo- 
nents of the four vector f£ are equal zero. If we study 
the matter composed of electrons, neutrons and protons, 
the zero-th component /" = fx is (see, e.g., Ref. [49j ) 



fx =V2G F 



E(A) 
n f q) , 

f=e,p,n 



lf } = (4l - 2Q (/) sin 2 9 W + 5 fe 5x„ e ) , (3.14) 

where n/ is the number density of background particles, 
is the third isospin component of the matter fcrmion 
/, is its electric charge, 6w is the Weinberg angle 
and Gf is the Fermi constant. 

Using Eqs. P^j) . (|3~3|) . and (153)) and we get that ma- 
trix (g a b) has the form, 



(9ab) 



f a cos 2 + fp sin 2 9 sin 28AV eS /2 
sin 26AV cS 12 fa sin 2 9 + fp cos 2 9 



(3.15) 



where AV c b = fp — f a is the difference between the ef- 
fective potentials of the flavor neutrinos interaction with 
background matter. With help of Eq. p. 141) we present 
A14ff in the form, 



AU cf 



V2G 



f x 



n e , for v e -> z/ M , T , 
0, for i/ MjT -> i/ Ti/t , 



(3.16) 
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for various oscillations channels. 

In the following we will discuss the low density matter 
limit, g a <C k, which is fulfilled for all realistic neutrino 
momenta and densities of background matter. Indeed, 
even for the background matter in the center of a neu- 
tron star where n n = 10 38 cm -3 , using Eq. p. 141) we get 
g a ~ 10 eV, which is much less than any reasonable neu- 
trino energy. With help of Eq. Q3.7[) we get that in this 
approximation w_/2 « $ — cos 29AV e s/2, where $ and 
Sm 2 are defined in Eq. (|2.2ip . 

The transition probability for the process vp — ¥ v a can 
be calculated on the basis of Eq. p. 131) as 

Pv,^u a (t) = K(z,t)\ 2 « P max sin 2 ( -^-t) , (3.17) 



where 

_ $ 2 sin 2 (2(9) 

max ~ (1 cos 29 - AV cS /2) 2 + $ 2 sin 2 (2(9) ' 

= ^($cos20- AV eB /2) 2 + $ 2 sin 2 (2(9), (3.18) 

-^osc 



are the maximal transition probability and the oscilla- 
tions length. 

One can see that Eqs. (|3.17|) and (|3.18|) reproduce the 
well known formula for the neutrino oscillations probabil- 
ity in the background matter (see Ref. [9(). If the back- 
ground density has the resonance value determined by, 
$ cos 29 — AV^f /2, the maximal transition probabil- 
ity reaches big values ~ 1. This resonance enhancement 
of neutrino oscillations in matter is known as the MSW 
effect @. 

Now let us discuss the modification of Eqs. (|3.17[) 
and (|3.18l) which include a hypothetical nonstandard 
neutrino interaction. One of the possibilities to account 
for this kind of interactions is to study the nondiagonal el- 
ement of the matrix (/\\/)- This interaction can produce 
the neutrino flavor conversion in presence of background 
matter. If still we discuss the nonmoving and unpolar- 
izcd matter, we define the additional nonzero element as 
/ = f°g 7^ 0. We can express the nonstandard interac- 
tion as / = e a p(fa + fp)/% although the exact form of 
the nonstandard interaction dependence on the densities 
of background fermions is still open. The experimental 
constraint on the e a p parameters reads \e a p\ i$ 0.4 [52j . 

Using the same technique as to get Eqs. (|3.17[) 
and (|3.18p we arrive to the modified maximal transition 
probability and the oscillations length, 



I 

P -= ^2B-^m^ + %^2B + fr £-^cos2,-AW2) 2 + ( $ sin2, + /) 2 . (3.19) 



Eq. ()3.19j) is exact and valid for arbitrary magnitude of 
the nonstandard interaction / in contrast to the pertur- 
bative formulas derived in Ref. 53]. 

As one can see in Eq. (13. 19)) that in the majority of 
cases the nonstandard neutrino interaction of the consid- 
ered type does not generate any additional resonances in 
neutrino oscillations. The small new interaction can just 
slightly change the shape of the transition probability. 

We can however notice that the new interaction pro- 
duces flavor oscillations even for massless neutrinos. In- 
deed, if we suggest that m a = (or $ = 0), we get that 
the parameters of transition probability formula, given 
in Eq. (|3.19j) . formally coincide with that in Eq. (|3.18|) . 
derived for the massive neutrinos, if replace $ sin 29 — > / 
there. However we cannot expect the appearance of the 
usual MSW resonance in this model since $ = 0. The 
amplification of neutrino oscillation can happen only if 
A14ff = 0. For example, it is the case for -H> v T oscil- 
lations [see Eq. fl3.16[) ]. 

Note that we have chosen the plane wave initial condi- 
tion corresponding to ultrarclativistic particles to study 
neutrino oscillations in background matter. It allowed 
us to exactly take into account the contribution of the 
field g in Eq. (|3.9[) . It is, however, possible to study the 



neutrino evolution with arbitrary initial condition in low 
density matter 34]. It was shown in Ref. [lij that the 
dynamics of neutrino oscillations is consistent with the 
results of Ref. [1] . 



IV. DIRAC NEUTRINOS IN AN EXTERNAL 
MAGNETIC FIELD 

In this section we apply the formalism developed in 
Sees. HT1 and HTT1 for the description of neutrino evolution 
in an external electromagnetic field [35]. In contrast to 
the previous sections we examine the situation when the 
helicity of a neutrino changes together with its flavor, 
i.e. we study so called neutrino spin flavor oscillations, 
^a ' L - We derive the new transition probability 
formulas which account for arbitrary magnetic moments 
matrix. 

Neutrinos are known to be uncharged particles. The 
constraint on the neutrino electric charge is at the level 
of 10~ 13 e [13] • Nevertheless there is a possibility for 
them to interact with an external electromagnetic field 
F^d = (E, B) via anomalous magnetic moments. The ex- 
perimental constraint on the neutrino magnetic moments 
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is ~ 10~ 10 /zb [H|,[54|, where /j,b is the Bohr magneton. 
Despite the smallness of magnetic moments, its inter- 
action with strong electromagnetic fields can produces 
sizeable effects (see, e.g., Sees. I VII and IVIll below) . 

The Lagrangian for the considered system of two flavor 
neutrinos v = (y a , i/p) is expressed in the following way: 

£ = ^ ^^l^d^vx 

A=a,/3 

- vximxx' + lMxxa^FvAvy, (4.1) 

\\>= a ,p \ z / 

where a^ v — (i/2)(j^ v — 7^7^). The magnetic moments 
matrix (M\\>) in Eq. (|4. 1[) is defined in the flavor eigen- 
states basis. In general case this matrix is independent 
from the mass matrix {mw), i.e. the diagonalization of 
the mass matrix does not necessarily imply the diagonal 
form of the magnetic moments matrix. 

To analyze the dynamics of the system (|4. 1[) we formu- 
late the initial condition problem [see Eqs. (|2.2p . (|2.6I) . 
and (|2.18[) ] and introduce the mass eigenstates ip a [see 
Eqs. (|2.3I) and (|2.4p ]. However, in contrast to the previ- 
ous sections we should choose the normalization spinor 
in Eq. (|2.18p in a specific form. 

When a neutrino with anomalous magnetic moment 
propagate in an external electromagnetic field its helicity 
changes. Therefore we can impose the additional con- 
straint on the initial spinor, 

M°> =■>?>, P ± =(l±£^>), (4.2) 

which means that one has neutrinos of the specific helic- 
ity initially. Here X = 7 s a is the Dirac matrix. If we 
act with the operator P T on the final state v a (r,t), we 
can study the appearance of the opposite helicity eigen- 
states among neutrinos of the flavor "a" , i.e this situa- 
tion corresponds to the neutrino spin flavor oscillations 
^a' L - For the sake of definiteness we choose the 
initial wave function corresponding to a left polarized 
neutrino and the final one to a right polarized particle. 

Now we express the Lagrangian (14.11) using the mass 
eigenstates ip a , which diagonalize the mass matrix, 

£ = Y ^aivy^dn - m a )ip a 
0=1,2 

~\ /WaCT^bF^, (4.3) 

afc=l,2 

where 

is the magnetic moment matrix presented in the mass 
eigenstates basis which, as we mentioned above, not nec- 
essarily to be diagonal. 



As in Sees. |TT] and IIIII we will study the situation of 
mass eigenstates neutrinos tf> a which are Dirac particles. 
It means that the magnetic moments matrix (fi a b) can 
have both diagonal and nondiagonal elements. The diag- 
onal elements of this matrix correspond to usual magnetic 
moments and the nondiagonal to the transition ones. The 
transition magnetic moments are responsible for the tran- 
sitions between left and right polarized particles of dif- 
ferent species. 

Let us assume that the magnetic field is constant, uni- 
form and directed along the z-axis, B = (0,0, B), and 
that the electric field vanishes, E = 0. In this case we 
write down the Pauli-Dirac equations for ip a , resulting 
from Eq. (|4.3|) . as follows: 

iipa =U a ip a +Vipb, a, 6 = 1,2, ffl^i, 

H a =(ap) + Pm a - /i a /3S 3 B, V = ~nPE 3 B, (4.5) 

where /i a = /i aa , and \i = fi a b = Hba are the elements of 
the matrix (// a &) (14.41) . 

We should notice that, as in Sec. IIIII the wave equa- 
tions (|4.5[) are coupled due to the presence of the inter- 
action V. Therefore we have to use a sort of the pertur- 
bative approach to account for this term, whereas anal- 
ogous diagonal magnetic interaction —ii a f3Y,zB will be 
taken into account exactly from the very beginning. 

We will study the propagation of neutrinos in 
the transverse magnetic field. Therefore it conve- 
nient to choose the initial momentum along the x- 
axis k = (k, 0,0) and the initial spinor as = 
(1/2)(1, —1, —1, 1). It is easy to see that the wave func- 
tion i/g\r) describes an ultrarelativistic particle prop- 
agating along the rr-axis with its spin directed opposite 
to the x-axis, i.e. a left polarized neutrino. The contri- 
bution of the longitudinal magnetic field to the dynam- 
ics of neutrino oscillations is suppressed by the factor 

m a /k<. 1 gnnm. 

The general solution of Eq. (|4.5[) can be presented as 
follows: 

x £ [ fl (0 (<)M (C) expH £(C) t) 
C=±i 

+ 6«)(t)«(«exp(+i£?(O t )]. (4.6) 

Our main goal is to determine the coefficients a^P and 
bcP consistent with both the initial condition (|2.18|) 
and l|4.2[) and the evolution equation (I4.5[) . As in Sec. IIIII 
these coefficients are in general functions of time. 

We have already mentioned that the helicity of a neu- 
tral particle with an anomalous magnetic moment is not 
conserved in an external magnetic field. Therefore to 
classify the states of massive neutrinos in Eq. (|4.6I) one 
has to use the operator [35j, [40J, |57| , 

n a = m a £ 3 + i 7 °7 5 (£ x p) 3 -y, a B, (4.7) 
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which commutes with the Hamiltonian T-L a in Eq. (I4.5[) 
and thus characterizes the spin direction with respect to 
the magnetic field. The quantum number Q ± 1 is the 
sign of the eigenvalue of the operator (I4.7[) . 

The energy levels in Eq. (|4.6[) have the form (3|1 H(| , 



eP 



pi 



c(C)2 



pi + pi- 



rn 



For our choice of the external magnetic field B = (0, 0, B) 
and the initial momentum k = (k, 0, 0), Eq. (|4.8|) reads 

2 

EP = JC a - C/iaS « k + ^ - Cfl a B, (4.9) 



for ultrarclativistic neutrinos with k 3> m a . Here JC a = 
\/k 2 + fflj is the kinetic energy of massive neutrinos. 

The exact form for the basis spinors uP and in 
Eq. (|4.6[) for arbitrary neutrino momentum is presented 
in Refs. [H, |57j- We reproduce the basis spinors for ul- 
trarelativistic neutrinos, 




(4.10) 



since we will be interested in the evolution of such par- 
ticles. In Eq. (|4.10|) we omit the index "a" since we ex- 
amine the case of k m a . The basis spinors u" and 
v~ correspond to the negative eigenvalue of the opera- 
tor (|4.7p (neutrino spin is directed oppositely to the mag- 
netic field), and u + and v + to the positive one (neutrino 
spin is parallel to the magnetic field). 

Using the general solution (|4.6p of the Pauli-Dirac 
equation (|4.5p containing the undetermined functions 



aP and bP and taking into account the orthonormality 



of the basis spinors (|4.10l) we get the system of ordinary 
differential equations for these functions aP 

iaP =exp(+iEPt)uP^V 



and bP: 



Eh 

C'=±i 



(C')JC') 



exp (-iE^'h) 



<7, 



exp i+iE^'k) 



ii>P =exp(-iEPt)vPv 



C'=±i 
.i,(C).,(0 



(C).,(C) 



exp (-iE^'h) 



exp (+iE { b c ' } t) 



(4.11) 



which should be supplied with the initial condition (|3.10p 
but with different initial wave function v\ = 
(1/2)(1, -1,-1,1) (see above). 



With help of the obvious identities (uf\V\u^) = Tl^B 
and (Ua\V\vf) = 0, which result from Eq. (|4. 10[) . one 
can cast Eq. (|4.1ip into the form 



it=TaffiBexp [i(£± - E±)t] 



(4.12) 



which is analogous to Eq. p. lip studied in Sec. Mil Note 
that the ordinary differential equations for the functions 



,(0 



and bP again decouple. 



On the basis of the results of Appendix [A] we are able 
to write down the solution of Eq. (|4.12l) as 



af(t) ^F ± af (0) + G ± af(0), 
o±(t) =F ± *af(0)-G ± *af(0), 



(4.13) 



where 



F ± 
G ± 



cos Q±t — i— — — sin Q±t 



2VL 



± 



exp (kj±t/2), 



± i— — sinf2-|-t exp (iw±i/2), 



and 



Ot = V(A^) 2 + (^±/2) 2 , w± = Ef 



Et 



(4.14) 



(4.15) 



The details of the derivation of Eqs. (14. 13l) - (|4. 15[) from 
Eqs. (|4.12p are also presented in Ref. [35| . 

Using Eq. (|4~^|) and Eqs. (|4TTT]) - (|4TT5]l and the identity 
[see Eq. (I4.10p ] we obtain the wave 
functions ip a , a = 1, 2, as, 

V'i (x, i) = exp PiE+t) (u + <8> u +t ) 

x [F+^iM + G+^faO)] 
+ exp (-LEft) ® w _t ) 
x [F-^i(x,0) + G-^ 2 (x,0)], 

ip%{x,t) = exp(-LE^i) (u + (giu +t ) 

x [F+> 2 (a;,0)-G + *Vi(^,0)] 



+ exp (— LE 2 t) (u ®u ') 
x [F~*ij; 2 (x,0)-G-*iJ 1 (x,0)} 



(4.16) 



which satisfy the chosen initial condition since G (0) = 
0, F ± (0) = 1 [see Eqs. (|4~14l) ] an d [(u+ ® u+t) + (it" ® 
«-t)]V» a (»,0) = ^a(»,0) [see Eq. gTU])]- 

To study the appearance of right polarized neutrinos 
of the type "a" we should act with the operator P + = 
(1 + Ei)/2 defined in Eq. (|4.2p on the final wave function 
v a (x,t), 

=^(1 + Si) 

x [cos6ip a (x,t) - sin6ip a (x,t)} , (4.17) 

where ip a (x,t) are shown in Eq. (|4.16l) . 
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With help of Eqs. (12^41 . (14T61) . and (14T71) we 

receive for the right polarized component of v a the ex- 
pression 



1 



sin o cos 



F-) 



+ cos 2 6»e- i,Clt (e iMlSt G + - e- iwBt G-) 
+ sin 2 ee-' lK2t {e^ 2Bt G + * - e~^ m G~*)} 

(4.18) 



where = (1/2)(1, 1, 1, 1) is the normalized 

spinor representing the right polarized final neutrino 
state. 



Finally, taking into account Eqs. (|4.14[) and (|4.15l) it 
is possible to express the wave function in Eq. (|4.18p in 
the form 



v^{x, t) =| sin 9 cos 9-^ 



sin(f2+i) exp (ifiBt) 



sin(J7_i) exp (—ifiBt) 



sin(i7+i) 9 
; cos 2 ( 



sin(f2_i) , 



cos(/i-Bi) 



x exp ( — icrt + ifcir)i^ 



(0)R 



(4.19) 



where (7 = (/Ci + /Ci)/2 and /2 = (/i! + /i2)/2- The mag- 
netic moments fi and /i a are defined in Eq. (I4.4[) . 

The transition probability for the process ffe — > v R can 
be directly obtained as the squared modulus of iy^(x,t) 
from Eq. (|4~T8|) or Eq. (|4~T9|) . that is P^l^r (i) = 

|z^(x, i)| 2 . Notice that the probability is a function of 
time alone with no dependence on spatial coordinates. 
This is of course obvious as we have taken the initial 
wave function as a plane wave and the the magnetic field 
spatially constant. 

Let us now apply the general results Eq. (|4.18j) or 
Eq. (|4.19l) to two special cases. We first consider the sit- 
uation where Ml, 2 3> /i, i.e. the case when the transition 
magnetic moment is small compared with the diagonal 
ones. Using Eqs. (|4.14p and (|4. 1 5[) we find that, in this 
case F ± « 1 and fl± ~ oj±/2, and Eq. (|4.18l) takes the 
form 

v%(x, t) wij sin 9 cos 9 [e~ iKlt sin mBt 

- e- iK2t sin fi 2 Bt] 

^ n ^B ( iV t sin A + t 
+ cos 26»— e~ lS+t — — - 
2 V A+ 



-iS_t 



sin A_i 



ifex (Q)R 



(4.20) 




vS(0 



FIG. 1. The schematic illustration of the system evolution 
in the case /i 2> Mi. 2- The horizontal lines of the figure cor- 
respond to various neutrino eigenstates at different moments 
of time (t — and t). The expressions next to arrows corre- 
spond to the appropriate factors in the formula (|4.21[) of the 
wave function. The arrows from f^(0) to ip\ (0) and ip2 (0), 
for example, indicate the vacuum mixing matrix transforma- 
tion at t = and the arrow from tp\ (0) to ^(t) the evo- 
lution of the mass eigenstates with the helicity change. The 
transitions ipa(Q) — ► can he described by the formula, 

^ 2 (() = Al^rV>2,i(0), where A L ->r = sin Q.t. This 

figure is taken from Ref. [35|. 



where 



S±=cr±^ J B, A ± =$±5fiB, 5/i = 



Mi - M2 



and the phase of vacuum oscillations <f> was defined in 
Eq. (|2~2ll . Eq. (|4~20l) was obtained in Ref. [H] using 
the perturbative methods. Assuming that /i <C /i a the 
perturbation theory was developed in that work. Now 
we rederive the same result as a particular case of the 
more general result. 

As another application of our general result we will 
study the situation, where the transition magnetic mo- 
ments are much larger than the diagonal ones, that is 
M ^ Mi. 2- In this case Eqs. (|4.14|) gives F + w F~ and 
G + w — G~ , and we receive from Eq. (|4.19p for the wave 
function i/ R the expression 



where 



v^{x, t) «i exp ( — icrt + \kx) cos(26*) 
x^ s in(QB«)^ 0)R , 



n B = v/(mS) 2 + * 2 - 



(4.21) 



(4.22) 



The transition probability for the process ijj — > v R is 
then given by 

^, S W=cos 2 (20)(^ sin 2 (^)- (4.23) 

The behavior of the system in this case is schematically 
illustrated in Fig [T] It should be noticed that the analog 
of Eq. (|4.23l) was obtained in Ref. [2(| where the au- 
thors studied the resonant spin flavor precession of Dirac 
and Majorana neutrinos in matter under the influence 
of an external magnetic field in frames of the quantum 
mechanical approach. 

Using Eq. (|4.19p one can describe spin flavor oscil- 
lations of Dirac neutrinos with arbitrary magnetic mo- 
ments matrix. It is the new result which was obtained 
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using the relativistic quantum mechanics approach. Nev- 
ertheless this result is consistent with the conventional 
quantum mechanical description of spin flavor oscilla- 
tions. We will demonstrate the consistency in Sec. [V] 
where more general case of neutrinos propagating in mat- 
ter and external magnetic field is studied. 

Spin flavor oscillations of Dirac neutrinos with arbi- 
trary initial condition, not necessarily corresponding to 
ultrarelativistic particles, were studied in Ref. |35j using 
the perturbative approach. To effectively apply the per- 
turbation theory one has to study the situation of small 
transition magnetic moment, \l -C fj, a - One could de- 
rive Eq. (|4.20[) using the results of Ref. [35| in the limit 
k > m a . 



V. DIRAC NEUTRINOS IN MATTER UNDER 
THE INFLUENCE OF A MAGNETIC FIELD 

In this section, using the relativistic quantum mechan- 
ics method, we study the general case of the mixed flavor 
neutrinos propagating in background matter and inter- 
acting with an external electromagnetic field [io| . We 
formulate the initial condition problem for neutrino spin 
flavor oscillations. Then we derive the effective Hamil- 
tonian which governs spin flavor oscillations and show 
the consistence of our approach to the usual quantum 
mechanics method. The corrections to the standard ef- 
fective Hamiltonian are also obtained. 

The Lagrangian for the system of two mixed flavor neu- 
trinos v = (u a ,vp) interacting with background matter 
and external electromagnetic field has the form, 



X=a,j3 XX'=a,/3 

+ ^M xx ,a^F^y x ,, (5.1) 



where the mass matrix (m\\>), matter interaction matrix 
if XX')) an d t ne magnetic moments matrix (M\x>) are de- 
fined in Sees. [TTl UTTl and IIVI respectively. 

In the following we will be interested in the stan- 
dard model neutrino interaction with matter which cor- 
responds to the diagonal matrix f£ x , = 5xx'f\- Moreover 
we will study the situation of nonmoving and unpolarized 
matter. In this case only zero-th component of the four 
vector /c is not equal to zero. The explicit form of this 
component f\ = /° for the background matter composed 
of electrons, protons, and neutrons is given in Eq. (|3.14|) . 

We choose the configuration of the electromagnetic 
field F^v = (E, B) in Eq. (|5.1[) in the same form as in 
Sec. IIVI Namely, we suppose that the electric field is ab- 
sent E = and magnetic field is constant and directed 
along the z-axis, B = (0, 0, B). 

To study the time evolution of flavor neutrinos we 
should supply the Lagrangian (|5.1I) with some initial con- 
dition. We choose the initial wave functions in the same 
form as in Sec. IIVI i.e. we suppose that u a (r, 0) = and 



v/3 (r, 0) = e lkx v^\ where the spinor uy corresponds to 
either left or right polarized neutrinos. The explicit form 
of Vp can be defined with help of the operators P± in 
Eq. (I4.2[) . It means that initially we have a beam of neu- 
trinos of the flavor "/3" with a specific polarization prop- 
agating along the x-axis. If we study the appearance of 
neutrinos of the flavor "a" of the opposite polarization, it 
will correspond to a typical situation of neutrino spin fla- 
vor oscillations in matter and transversal magnetic field, 
l,r , . R l 

Then we introduce the mass eigenstates ip a using 
Eqs. (|2.3p and (|2 .4[) to diagonalize the mass matrix 
(mxX') m Eq. (|5.ip . These mass eigenstates are again 
supposed to be Dirac particles. Now the Lagrangian (|5.1|) 
expressed via the new mass eigenstates has the form, 



£ = E V'aOV^ ~ «1a)-0a 

L,2 

E & {dubll + Mv^fA fa, (5.2) 



o=l,2 



ab=l,2 v ' 

where (g^ b ) and (p a b) are the matrix of neutrino inter- 
action with matter and the neutrino magnetic moments 
matrix expressed in the mass eigenstates basis, which are 
defined in Eqs. (|3.3p and (|4.4|) . We remind that in case 
of nonmoving and unpolarized matter the matrix (g^ b ) 
has only zero-th component (I3.5|) . 

On the basis of the mass eigenstates Lagrangian (|5.2[) , 
we can derive the corresponding wave equations which 
have the following form: 

i^o =^a^a + Vipb> a = 1,2, a^b, 
H a =(ap) + f3m a - ^ a /3^ 3 B + g a (l - 7 5 )/2, 
V = -/i/3S 3 S + .g(l- 7 5 )/2. (5.3) 

Note that we cannot directly solve the wave equa- 
tions (|5.3p because of the nondiagonal interaction V 
which mixes different mass eigenstates (see also Sees. IIIII 
and llVp . Nevertheless we can point out an exact solution 
of the wave equation iip a = H a ipa, for a single mass eigen- 
state ipa, that exactly accounts for the influence of the 
external fields g a and ^i a B. The contribution of the mix- 
ing potential V can be then taken into account using the 
perturbation theory, with all the terms in the expansion 
series being accounted for exactly. 

We look for the solution of Eq. (5.3) in the following 
form EH: 



d 3 p 



= ipr 



(2tt) 3 / 2 
[aP(t)uPexp(-iEPt) 



C=±i 



(5.4) 



where the energy levels, which were found in Ref. [37j |. 
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have the form, where R 2 = \fp 2 M^ + (fJ, a B) 2 m a . 



The basis spinors in Eq. (|5.4|) can be found in the limit 



M a -\J (HaB) 2 + ejf/J, (5.5) of the gmall neutrin0 mass ^3, 
I 



t (0 



,(0 



2^2M a {M a - Cffa/2) 



2^/2Ma(M a +(ga/2) 



( H a B + CMa- 9a/T 
fJL a B - (M a + 9a/2 
H a B - (M a + 

\fJ, a B + (Ma - 9a/2 

M a -CWB-g a /2] 
M a + {[fi a B + g a /2] 
-M a -{[fi a B + g a /2] 
-M a + Chi a B - g a /2] 



(5.6) 



It should be noted that the discrete quantum number 
C = ±1 in Eqs. (|5.4[) - (I5.6[) does not correspond to the 
helicity quantum states. 

Now our goal is to find the time dependent coeffi- 
cients do (t) and bu(f). On the basis of the general 
solution (|5.4p of the wave equation (|5.3[) we obtain the 
ordinary differential equations for these functions which 
formally coincide with Eq. (|3.9p . However the mixing 
potential V is now defined in Eq. (|5.3|) . To obtain the 
modified Eq. (|3.9I) we again use the orthonormality of the 
basis spinors (|5.6I) . The initial condition for the functions 

a a C \t) and b a C \t) also coincide with Eq. (pTTO]) . with vf ] 
corresponding to a definite helicity spinor. 



Taking into account the fact that (u^ \V\vjf ') = 0, 

we get that the equations for da (t) and b£\t) decouple, 
i.e. the interaction V does not mix positive and negative 
energy eigenstates. In the following we will consider the 

evolution of only da (t) since the dynamics of ba(t) is 
studied analogously. 

Let us rewrite the modified Eq. (|3.9p in the more 
conventional effective Hamiltonian form. For this 
purpose we introduce the "wave function" ^ T = 
(a^, o^, , ). Directly from the modified Eq. (|3 .9[) 
for the functions a a ^(t) we derive the equation for 



d*' 

i^f - = JZV, H' = 
At 



( 


h- 







H _(P-t\ 


h-e-' 1 "-* 

































/ 



(5.7) 



where 



h T =(u^\V\uf 



1 



8y/M a Mb(M a ± 9a/2)(M b ± 5 ft/2) 

x [2^B{g a ^ b B + .g 6 Aia-B) ± A^iB(n a BM b + ^ b BM a ) ± 2g(g a M b + g b M a ) + 4gM a M b + gg a g b ], a^b, 

1 



ff T =(«?|V|4) = (uf\V\uf) 



±1 



V-MiW-Mi ± gi/2){M 2 T 52/2) 
x [2^iB(g 1 fj, 2 B + g2^iB) qp 4(iB(fixBM2 - H2BM1) T 2g(giM2 - gzMi) - ^gMiM.2 + 951^2] 



(5.8) 



as well as w T = Ef - Ef + (51 - 0j)/2 and ft T = £f - £^ + (pi - ,g 2 )/2. 

Instead of it is more convenient to use the transformed "wave function" "J defined by 



=UV, W = diag|e i 



(n+w_)t/2 i(Q—U)-)t/2 — i(f2-w + )t/2 -i(fi+w + )t/2 



}■ 



(5.9) 



where 17 = (f2_ — f2_|_)/2, to exclude the explicit time dependence of the effective Hamiltonian H'. Using the property 
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lu + + cj_ = Sl + + f2_, we arrive to the new Schrodinger equation for the "wave function" 'J, 



d* 



H = WH'U - tLVU 



/{SI 


+ w- 


_)/2 h. 







h- 


(0-w_)/2 




\ 







-(SI-oj 


H_ 










-(fi + w+)/2 



(5.10) 



Despite initially we used perturbation theory to account for the influence of the potential V on the dynamics of the 
system (|5.3|) , the contribution of this potential is taken into account exactly in Eq. (|5.10|) . It means that our method 
allows one to sum up all terms in the perturbation series. 

As we mentioned above, the quantum number £ does not correspond to a definite helicity eigenstate. Thus the 
initial condition, which we should add to Eq. (|5.10l) . has to be derived from Eqs. (|3. 10[) and (|5.6[) and also depend 
on the neutrino oscillations channel. For example, if we discuss —> neutrino oscillations, the proper initial 
condition for the "wave function" \l/(0) = \&o is 



= I —sm0\l , — cos0\j ,sm0\l ,cosP 



2Mi 



2M, 



2Mi 



2M, 



(5.11) 



Suppose that one has found the solution of the system (I5.10[) and (|5.11[) as ^(t) — (ipi,ip2,ip3,ip4)- Then the 
transition probability for v\ — > oscillations channel can be found as 



v n\ 1 J Mi^cosfl 

K^ufW =2] 



{ 4Ml 

H2B sin 8 



/M 7 



^Mx+ gi /2 y/M!-gi/2 



VM2+92/2 ^JM.2 - 52/2 



(5.12) 



To obtain Eq. (|5.12|) for simplicity we use the fact that 
initially we have rather broad (in space) wave packet, 

corresponding to the initial condition vp(v, 0) = e lkx Vp 
[see Eq. pT8]) ]. 

Now we demonstrate the consistency of the results of 
relativistic quantum mechanics approach to the descrip- 
tion of neutrino spin flavor oscillations (see Eqs. (|5.10[1 - 
(I5.12p . which look completely new) with the standard 
quantum mechanical method developed in Ref. [20|. We 
remind that the following effective Hamiltimian: 



H'qm — 



* + 9i 


9 




9 


-$ + 92 


-HB 


-HiB 




$ 


-fiB 


-H2B 







(5.13) 



was proposed in Ref. [20| to describe the evolution of 
neutrino mass eigenstates in matter under the influence 
of an external magnetic field. 

The effective Hamiltonian H'q M acts in the space with 
the basis composed of helicity eigenstates of massive neu- 
trinos. As we mentioned above, the helicity operator 
(Sp)/|p| does not commute with the Hamiltonian H a in 
Eq. (|5.3p . Therefore the choice of the helicity eigenstates 
as the basis functions is justified only in the relatively 



weak external mag netic field case (see the detailed dis- 
cussion in Ref. [40]) or in case of the small diagonal mag- 
netic moments [39]. In our approach we use the basis 
spinors (|5 .6[) which are the eigenfunctions of the Hamil- 
tonian Ha and exactly take into account matter density 
and magnetic field strength. Thus these spinors are more 
appropriate basis functions for the description of spin fla- 
vor oscillations. 

We have found that in frames of the relativistic quan- 
tum mechanics approach the dynamics of the neutrino 
system can be described by the Schrodinger like equation 
with the effective Hamiltonian (|5.10[) . Let us decompose 
the energy levels (I5.5[) supposing that neutrinos are ul- 
trarelativistic particles, 



eP = k + 



9a 

2 



+ c 



2k ^8k 2 M a 



(5.14) 



In Eq. f|5 . 14[) we keep the term ~ m a /k 2 to examine 
the corrections to the conventional quantum mechanical 
approach. 

Performing the similarity transformation of the effec- 
tive Hamiltonian H in Eq. (|5.10[) and using the orthogo- 
nal matrix 1Z (1Z T 1Z = I) of the following form: 
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n 



( ^M 1 +g l /2 

o 

\ 





y/M 2 + g2/2 
V2~Mz~ 



- <? 2 /2 



Mi5 



y/2Mi(Mi+ gi /2) 


y/2Mi{Mi- gi /2) 






M 2 S 



V27W 2 (^2 + .g 2 /2) 


/i 2 S 



y27R; y/2M 2 {M 2 -g2/2) / 

we can see that the Hamiltonian H transforms to 1Z T H1Z ss Hqm + SH, where 



H, 



QM 



/$ + 3.gi/4- g 2 /4 g - f nB -fj,B 

g -$ + 3.g 2 /4-.9i/4 —fj,B -\i 2 B 

-HiB -liB $ - (.gi + ,g 2 )/4 

V -liB -» 2 B (,9i +.g 2 )/4 



and 



16/c 2 



2 fli 



,2 5I 2 gf 2 g% 



SH = ^dmg ( -™i^|.-™ 2 ^ 1 %^,m 2 ^ 



(5.15) 



(5.16) 



(5.17) 



is the correction to the standard effective Hamiltonian. 
It should be noted that the transformation matrix 1Z in 
Eq. (I5.15[) depends on the magnetic field strength and 
the matter density. 

The effective Hamiltonian Hqm is equivalent to Hq M 
in Eq. (|5.13j) since H QM = H' QM - tx(H' QM )/A ■ /, where 
I is the 4x4 unit matrix. It is known that the unit 
matrix does not change the particles dynamics. Thus the 
relativistic quantum mechanics approach is equivalent to 
the standard approach developed in Ref. [2(| • 

Now let us discuss the correction SH [see Eq. (|5.17l) ] 
to the quantum mechanical method. This correction re- 
sults from the fact that we use the correct energy levels 
for a neutrino moving in dense matter and strong mag- 
netic field. Note that in Eq. (15.171) we keep only the 
diagonal corrections ~ rn 2 /k 2 to the effective Hamilto- 
nian (|5.16p . If we slightly change non-diagonal elements 
of the Hamiltonian, it will result in the small changes 
of the transition probability. However, if we add some 
small quantity to diagonal elements, it can produce the 
resonance enhancement of neutrino oscillations. 

We should remind that the expressions for the basis 
spinors (|5.6p were obtained in the approximation of neu- 
trinos with small masses, whereas in Eq. (|5.14[) we ex- 
pand the energy up to ~ m 2 /k 2 terms. If we take into 
account ~ m 2 /k 2 corrections to the basis spinors (|5.6I) . 
we can expect that some non-diagonal entries in the ef- 
fective Hamiltonian H (|5.10[) will also obtain ~ mt/k 2 
corrections: h± — > h±+Sh± and H± — > H±+SH±. How- 
ever, using the explicit form of the effective Hamiltonian 
H (|5.10l) and the matrix 7Z (j5 . 15[) we get that these ad- 
ditional contributions are washed out in diagonal entries 
in Eq. (|5.17[) . We analyze the validity of the approxima- 
tions made in the derivation of the correction (|5.17l) in 



Appendix [B] 



VI. SPIN FLAVOR OSCILLATIONS OF DIRAC 
NEUTRINOS IN THE MAGNETIZED 
ENVELOPE OF A SUPERNOVA 

In this section we study the application of the general 
formalism for the description of neutrino spin flavor os- 
cillations, developed in Sec. [V] to the situation of neutri- 
nos propagating in the expanding envelope formed after 
a supernova explosion [391 ] . We find an exact solution 
of the Schrodinger equation with the effective Hamilto- 
nian (|5.16p for the background matter profile present in 
an expanding envelope and in the supernova magnetic 
field. We also analyze the possibility of enhancement of 
neutrino oscillations. 

To describe the dynamics of neutrino spin flavor oscil- 
lations one has to solve the evolution equation with the 
Hamiltonian (|5.16[) . This problem, in its turn, requires 
to solve a secular equation which is the fourth-order al- 
gebraic equation in order to find the eigenvalues of the 
effective Hamiltonian. Although one can express the so- 
lution to such an equation in radicals, its actual form 
appears to be rather cumbersome for arbitrary parame- 
ters. 

If we, however, consider the case of a neutrino prop- 
agating in the electrically neutral isoscalar matter, i.e. 
n e = n p and n p = n n , a reasonable solution is pos- 
sible to find. We will demonstrate later that it corre- 
sponds to a realistic physical situation. As one can in- 
fer from Eq. f|3 . 14[) for the case of the v\ — > oscilla- 
tions channel, in a medium with this property one has 
the effective potentials f a = = = — Gf?i/\/2 an d 



1G 



-.92 = go, where 
G F n/V2. 



fp = fe — Ve — G F n/y/2, where n 
Using Eq. (|3.15|) we obtain that g\ - 
g Q = -Ucos26», g = V sin 20, and V = 

Let us point out that background matter with these 
properties may well exist in some astrophysical environ- 
ments. The matter profile of presupernovae is poorly 
known, and a variety of presupernova models with differ- 
ent profiles exist in the literature (see, e.g., Ref. .58]). 
Nevertheless, electrically neutral isoscalar matter may 
well exist in the inner parts of presupernovae consisting 
of elements heavier than hydrogen. Indeed, for example, 
the model W02Z in Ref. [58[ predicts that in a 15Mq 
presupernova one has Y e — n e /(n p + n n ) = 0.5 in the 
O+Ne+Mg layer, between the Si+O and He layers, in 
the radius range (0.007-0.2)i? Q . 

We also discuss the model of neutrino magnetic mo- 
ments in which the nondiagonal elements of the magnetic 
moments matrix (|4.4[) are much bigger than the diagonal 
magnetic moments. Such a magnetic moments matrix 
was previously discussed in our works [35|,[33,[39| ( see a ^ so 
Sec. ITV)) . Note that in case of negligible diagonal mag- 
netic moments the helicity operator (|4.2[) commutes with 
the Hamiltonian H a in Eq. (|5.3j) and hence the effective 
Hamiltonian (|5.16j) acts in the helicity eigenstates basis. 
In other words the effective Hamiltonian (|5.13|) proposed 
in the standard quantum mechanical approach [20] is jus- 
tified in our case. 

For neutrinos having such magnetic moments and 
propagating in isoscalar matter the effective Hamilto- 
nian (|5.16[) is replaced by 



H, 



QM 



+ go 
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V -M-B 
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-fi,B 


+ go) 


-\iB 
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(6.1) 



We now look for the stationary solutions of the 
Schrodinger equation with this Hamiltonian. After a 
straightforward calculation one finds 



C=±i 

(v c <8> V/) exp (i£ c t)] * , 



where we have denoted 



£± =-^j2V 2 +4(^B) 2 + 4$ 2 - 4$U cos 29 ± 2VR, 



R = y/(V- 2$ cos 26») 2 + 4(fiB) 2 . 



(6.3) 



The vectors U± and V± are the eigenvectors correspond- 
ing to the energy eigenvalues £± and —£±, respectively. 



They are given by (£ = ±) 



z< 

sin20(£ c - $) 

-fiB sin 29 
-VBZ C /(£ C + 4>) 

- sin 29(£ c - $) 



c 



where 



V + (R 



iv c 2 



1 + 



fiBZ c /(£ c + $) 
-fiB sin 20 



- £(; COS 20 

(/x£) 2 



(6.4) 



$)2 

\2 



+ sin 2 (29) [IjiBf + (£ C $) 2 ] . 



(6.5) 



It should be noted that Eq. (I6.2p is a general solution 
of the evolution equation with the effective Hamilto- 
nian (16.11) satisfying the initial condition ^(0) = Vto- 

Note that we received the solution (|6.2[) - (|6.5[) under 
some assumptions on the external fields such as isoscalar 
matter with constant density and constant magnetic 
field. In Sec. [V] we showed that our method is equiva- 
lent to the quantum mechanical description of neutrino 
oscillations [2(| which can be used for a more general case 
of coordinate dependent external fields. Nevertheless the 
assumption of constant matter density and magnetic field 
is quite realistic for certain astrophysical environments 
like a shock wave propagating inside an expanding enve- 
lope after a supernova explosion. 

Consistently with Eqs. (|2~2"1) - (ET4]) . and (|2TT5)) with an 
initial spinor corresponding to a left polarized neu- 
trino (14.21) . we take the initial wave function ^(0) = in 
Eq. dnH) as *J = (^l^l^r ^R) = ( S i n 6), cos 6», 0,0). 
Using Eqs. (|6.2|) - (|6.5[) one finds the components of the 
quantum mechanical wave function corresponding to the 
right polarized neutrinos to be of the form 



Zl 



sin 2 (20)(£+ - <S>)e- i£+t 



sin 9 sin 29Z, 



-iS+t 



£+ + <$> 



+ {+^ -}, 



= 4 (sin0 



sin 2 (26») (£+ - $)e 



iS+t 
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cos 9 sin 29 Z + 



AE+t 



-i£+t£+ 



+ {+^ -}, 



(6.6) 



where the {+ — > — } stand for the terms similar to 
the terms preceding each of them but with all quan- 
tities with a subscript + replaced with corresponding 
quantities with a subscript — . The wave function of 
the right-handed neutrino of the flavor "a", v^, can 
be written with help of Eqs. (|2T2" ]) -([2"^ )> . and (|61)f as 
i^(t) = cos^VfW - sin6»V^ (t). 

The probability for the transition ih — > is obtained 
as the square of the quantum mechanical wave function 
i/;?. One obtains 



where (£ = ±) 



= [C+ cos(£+t) + C- cos(£_^)] 2 
+ [S+ sin(£+t) + S- sin(£_i)] 2 , 



(6.7) 



sin 2 (2tf)(£ c -<I>) , 



cos 29 



8 C + $ 



sin 2 (260(£ c -$) 



(6.8) 



As a consistency check, one easily finds from Eq. 
that C + + C_ = as required for assuring P(0) = 0. 

In the following we will limit our considerations to the 
case £+ ~ corresponding to the situations where 
the effect of the interactions of neutrinos with matter 
(V) is small compared with that of the magnetic interac- 
tions (fJ.B) or the vacuum contribution (<£>) or both [see 



Eq. (16. 3[> ]. Note that in this case one can analyze the ex- 
act oscillation probability (I6.7[) analytically, which would 
be practically impossible in more general situations. 

In the case ~ one can present the transition 
probability in Eq. (|6.7p in the following form: 



P{t) = P Q (t) + P c (t) cos(2m) + P a (t) sin(2m), (6.9) 



where 



p o(t) =2 [ s + + S-+ 2S+S- cos(25fii) 
™4C* + C_sin 2 ((5nt)], 



S 2 _) cos{2dnt) + 2S+S- 
4C+C-sm 2 {5nt)], 



Ps{t) =-{S 2 + - S 2 _)sm{2Snt), 



(6.10) 



and 



n = 



£+ + £- 



sn 



£+-£- 



3.11) 



As one can infer from these expressions, the transition 
probability P(t) is a rapidly oscillating function, with the 
frequency f2, enveloped from up and down by the slowly 
varying functions P u ,d — Po ± \/P 2 + P 2 , respectively. 

The behavior of the transition probability for various 
matter densities p and the values of \xB and for a fixed 
neutrino energy of E — 10 MeV and squared mass differ- 
ence of Sm 2 = 8 x 10~ 5 eV 2 is illustrated in Figs. [2]|4j 

As these plots show, at low matter densities the en- 
velope functions give, at each propagation distance, the 
range of the possible values of the oscillation probability. 
At greater matter densities, where the probability oscil- 
lates less intensively, the envelope functions are not that 
useful in analyzing the physical situation. 

One can find the maximum value of the upper envelope 
function, which is also the upper bound for the transition 
probability, given as 



(S+-S-) 2 , iiB<B', 

^i max) = <! C+C-jSl-Sl) 2 (6.12) 

C+C_(5 2 +5 2 ) + (C+C_) 2 + (5+5_) 2 ' 1 



where the value B' is the solution of the transcendent 
algebraic equation, C + C_ = S + S-. The corresponding 
maximum values of the averaged transition probability 
Pq{x) are given by 

pirn**) = 1 [(s+5 _)2 _ 4C+C-], (6.13) 

for arbitrary values of B. The values of these maxima de- 
pend on the size of the quantity fiB. These dependencies 



are plotted in Figs. (Ud)-|3|d). In the case of rapid oscil- 
lations the physically relevant quantities, rather than the 
maxima, are the averaged values of the transition prob- 
ability, which are also plotted in these figures. 

As Figs. [2jd)-|4|d) show, the interplay of the matter 
effect and the magnetic interaction can lead, for a given 
magnetic moment fi, to an enhanced spin flavor transition 
if the magnetic field B has a suitable strength relative 
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FIG. 2. (color online) (a)-(c) The transition probability versus the distance passed by a neutrino beam in matter with the 
density p = lOg/cc; (a) pB = 5 x 10" 13 eV, (b) pB = 1.1 x 10" 12 eV, (c) pB = 5 x 10" 12 eV. We take that E v = lOMeV, 
Sm 2 = 8 x 10 -5 eV 2 and 8 = 0.6, which is quite close to the solar neutrinos' oscillations parameters. The black line is the 
function P(x), the blue and green lines are the envelope functions P u ,d(x), and the red line is the averaged transition probability 
Po(x). (d) The dependence of the maximal values of the functions P{x) and Po(%), blue and red lines, respectively, on the 
magnetic energy pB for the given density. This figure is taken from Ref. [3Sj | - 



to the density of matter p. In our numerical examples 
this occurs at pB mllx = 1.1 x 10~ 12 eV for p = lOg/cc, 
at pB niax = 6.6 x 10~ 13 cV for p = 50g/cc, and at 
M-B max = 8 x 10~ 13 eV for p = lOOg/cc. For these values 
of pB both the maxima and the average of the transi- 
tion probability become considerably larger than for any 
other values of pB. Figs.[5{b)-|3|b) correspond to the sit- 
uation of maximal enhancement, whereas Figs. [Ha)-|Ua) 
and Figs. [UJc)-[4jc) illustrate the situation above and be- 
low the optimal strength B max of the magnetic field. 

It is noteworthy that the enhanced transition proba- 
bility is achieved towards the lower end of the pB re- 
gion where substantial transitions all occur, that is, at 
relatively moderate magnetic fields. At larger values 
of pB the maximum of the transition probability ap- 
proaches towards cos 2 (2#). Indeed, if pB ^> m&x(&,V), 
the transition probability can be written in the form 
P(t) = cos 2 (29) sin 2 (pBt) [see Eq. ([£25]) ]. It was found 
in Ref. [HH that neutrino spin flavor oscillations can be 
enhanced in a very strong magnetic field, with the tran- 



sition probability being practically equal to unity. This 
phenomenon can be realized only for Dirac neutrinos with 
small nondiagonal magnetic moments and small mixing 
angle. As we can see from Figs. [H[d)-|3Id) the situation 
is completely different for big nondiagonal magnetic mo- 
ments. 

One should notice that for long propagation distances 
consisting of several oscillation periods of the envelope 
functions, the enhancement effect would diminish consid- 
erably because of averaging. In the numerical examples 
presented in Figs. [2j4] the period of the envelope function 
is of the order of 10 3 — 10 4 km, which is a typical size of 
a shock wave with the matter densities we have used in 
the plots (see, e.g., Ref. [6fjj). Thus the enhanced spin 
flavor transition could take place when neutrinos traverse 
a shock wave. 

Let us recall that the above analysis was made by as- 
suming neutrinos to be Dirac particles. We will see below 
(see Sec. [X]) that the corresponding results are quite dif- 
ferent in the case of Majorana neutrinos. 
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FIG. 3. (color online) The same as in Fig. [2]for the density p 
liB = 5x 1(T 12 eV. After Ref. HH. 

VII. SPIN FLAVOR OSCILLATIONS BETWEEN 
ELECTRON AND STERILE ASTROPHYSICAL 
NEUTRINOS 

In this section we continue with the studies of oscilla- 
tions of supernova neutrinos on the basis of the effective 
Hamiltonian derived in Sec. [Vj In particular, we discuss 
the possibility of spin flavor oscillations between right po- 
larized electron neutrino and additional sterile neutrino 
in an expanding envelope formed after a supernova ex- 
plosion under the influence of a strong magnetic field. It 
is shown that the resonance enhancement of neutrino os- 
cillations is possible if we take into account the correction 
to the effective Hamiltonian (|5.17[) . 

As we mentioned in Sec. IVI1 in general case the solu- 
tion of the Schrodinger equation based on the effective 
Hamiltonians (I5.16|) and (|5.17[) is quite cumbersome and 
should be analyzed numerically. To study spin flavor os- 
cillations analytically we consider, just for simplicity, the 
situation when the vacuum mixing angle between differ- 
ent neutrino eigenstates is small, 9 <C 1. If 6 -C 1, the 
matter interaction term and the the magnetic moments 
coincide in both flavor and mass eigenstates bases. 

A resonance of neutrino oscillations can appear if the 
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difference between two diagonal elements in the effective 
Hamiltonian is small [6lj . In table HI where we take into 
account the contributions of both Eqs. (|5.16p and (|5.17l) . 
we list the resonance conditions for various oscillations 
channels. It can be seen from this table that the correc- 
tions (|5.17p to the effective Hamiltonian become impor- 
tant when we discuss vf ^ T O ,R oscillations channel, 
where v s is the additional sterile neutrino. 

As an example, we can study the case of vf f-> v\ 
oscillations. Putting a = s and /3 = e in table U and 
using Eq. (|3.14p . we obtain that f s = 0, f e = V2Gp(n e — 
n n /2) and M 2 = M e = From table U 

we obtain the resonance condition for this oscillations 
channel as, 

where we discuss the case of small diagonal magnetic mo- 
ment of an electron neutrino: [i Vc B <C f e and E v = k is 
the neutrino energy. In Appendix [B] we discuss other 
small factors which can also contribute to the resonance 
condition (17. ip . 

As an example we consider -o- v\ oscillations in 
an expanding envelope of a supernova explosion. It is 
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FIG. 4. (color online) The same as in Fig. [2] for the density p = lOOg/cc; (a) pB = 4 x KT 13 eV, (b) fiB = 8 x 10~ 13 eV, (c) 
fiB = 5x 1CT 12 eV. After Ref. 

TABLE I. The resonance conditions for the -O Vg oscillations channels which take into account the Hamiltonians (|5.16|) 
and ((5T7T) . 



No. 



Oscillations channel 



Resonance condition 



Va -S-S- v n 
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2 32fc 2 



ml fa 



+ 



"2 



//3 - fa 



32k 2 
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32k 2 

m\fl 
Mj 



Mj 

mjfl 
Mj 

M\ 



M\ 



known (see, e.g., Refs. [(32|,|63| for the most detailed anal- 
ysis) that right polarized electron neutrinos can be cre- 
ated during the explosion of a core-collapse supernova. 
Indeed, if a neutrino is a Dirac particle, the spin flip 
can occur in the following reaction: + (e~ , p, N) — > 
vr + {e >P, N), with electrons e~ , protons p, and nuclei 
N in the dense matter of a forming neutron star. This 



spin flip is due to the interaction of the neutrino diagonal 
magnetic moment with charged particles. Hence left po- 
larized neutrinos are converted into right polarized ones 
with energy in the range (100 - 200) MeV fo4} . 

Besides the fact that the creation of right polarized 
neutrinos can provide additional supernova cooling [6o| . 
these particles can be observed in a terrestrial detec- 
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tor. When right polarized neutrinos propagate from the 
supernova explosion site towards the Earth, interacting 
with the galactic magnetic field, their helicity can change 
back and they become left polarized neutrinos. Although 
the flux of this kind of neutrinos is smaller than that of 
generic left polarized neutrinos [63|, [66j], potentially we 
can detect these particles. The analysis of the neutrino 
spin precession in the galactic magnetic field is given in 
Ref. m- 

Suppose that the flux of right polarized electron neu- 
trinos is crossing an expanding envelope of a supernova. 
A shock wave can be formed in the envelope [67] . Ap- 
proximately 1 s after the core collapse, the matter den- 
sity in the shock wave region L ~ 10 8 cm can be up to 
10 6 g/cm 3 . We can also suppose that the matter density 
is approximately constant inside the shock wave (see also 
Sec. EQ- 

For the electroneutral matter Eq. (|7.ip reads 
8m 2 «5.0 x l(T 17 eV 2 x (3Y e - 1) 




where Y e = n e /(n e + n n ) is the electrons fraction. From 
Eq. (|7.2I) we can see that for the matter with Y e > 1/3 
and p ~ 10 6 g/cm 3 and an electron neutrino with E v ~ 
100 MeV and m„ c ~ 1 eV the mass squared difference 
should be 5m 2 ~ 10" 17 eV 2 . 

Note that the possibility of existence of sterile neutri- 
nos closely degenerate in mass with active neutrinos was 
recently discussed. In Ref. [68] it was examined how the 
fluxes of supernova neutrinos can be altered is presence 
of the almost degenerate neutrinos. The effect of spin fla- 
vor oscillations between active and sterile neutrinos with 
small dm 2 on the solar neutrino fluxes was discussed in 
Ref. [69(. The implications of the CP phases of the mix- 
ing between active and sterile neutrinos in the scenario 
with small Sm 2 to the effective mass of an electron neu- 
trinos were considered in Ref. [7(j. In Ref. 71] it was 
examined the possibility to experimentally confirm, e.g., 
in the IceCube detector, the existence of almost degener- 
ate in masses sterile neutrinos emitted with high energies 
from extragalactic sources. The range of 5m 2 studied 
in Ref. is 10" 16 eV 2 < Sm 2 < 10~ 12 eV 2 which is 
quite close to the estimates obtained from Eq. (I7.2[) . In 
Ref. [7l| the sterile neutrinos with even smaller dm 2 were 
studied: 10" 19 eV 2 < Sm 2 < 10- 12 eV 2 . 

Besides the fulfillment of the resonance condition (|7.ip , 
to have the significant o v\ transitions rate the 
strength of the magnetic field and distance traveled by 
the neutrinos beam, which we take to be equal to the 
shock wave size L, should satisfy the condition \xBL « 
7r/2. We can rewrite this condition as, 

For the transition magnetic moment /i = 3 x 10~ 12 /zb [72[ 



and L ~ 10 s cm (see above), we get that B ~ 10 7 G. Sup- 
posing that the magnetic field of a neutron star depends 
on the distance as B(r) = BQ(R/r) 3 , where R = 10km 
is the typical protoneutron star radius and Bq = 10 13 G 
is the magnetic field on the surface of a protoneutron 
star, we get that at r = 10 8 cm the magnetic field reaches 
10 7 G, which is consistent with the estimates of Eq. (17.31) . 
Note that magnetic fields in a supernova explosion can 
be even higher than 10 13 G and can reach the values of 
- 10 16 G [73|. 

VIII. MAJORANA NEUTRINOS IN VACUUM 

In this section we apply the formalism, developed in 
Sees, imrvl for the description of the Dirac neutrinos evo- 
lution in various external fields, for the studies of oscilla- 
tions of Majorana neutrinos. Here we study the evolution 
of mixed Majorana neutrinos in vacuum. We solve the 
initial condition problem for the two component neutrino 
wave functions and find the transition probability. How- 
ever, besides flavor oscillations, for Majorana particles 
one can consider transitions between neutrinos and an- 
tineutrinos. We study this process also on the basis of 
relativistic quantum mechanics. 

In general case the system of flavor neutrinos can be 
described with help of the appropriate number of left 
and right handed spinors (44}. We can however sug- 
gest that only left handed chirality components v\ = 
(1/2)(1 — 7 5 )^a are present in our system. It is the 
case, for example, for neutrinos in frames of the stan- 
dard model where the right-handed neutrinos are sterile. 
The general situation when one has both left and right 
handed particles will be studied in Sec. IXII 

The general mass matrix, involving left-handed flavor 
neutrinos v\ , can be diagonalized with help of the matrix 
transformation [44j . 

a 

where A = a,/3 is the flavor index and r] a , a = 1,2, 
corresponds to a Majorana particle with a definite mass 
m a . In the simplest case the mixing of the flavor states 
arises purely from Majorana mass terms between the left- 
handed neutrinos, and then the mixing matrix U\ a is a 
2x2 and unitary matrix, i.e., a = 1, 2 and, assuming no 
CP violation, it can be parameterized in the same way 
as in Eq. ([2^1 . 

We study the evolution of this system with the follow- 
ing initial condition [see also Eq. (|2.18p ]: 

i£(r,0)=0, ^(r,0) = ^ O) e ikr , (8.2) 

where k = (0, 0, k) is the initial momentum and v\ — 
(0, 1). The initial state is thus a left-handed neutrino of 
flavor "/?" propagating along the z-axis to the positive 
direction. 
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As both the left-handed state v\ and Majorana state 
rj a have two degrees of freedom, we will describe them in 
the following by using two-component Weyl spinors. The 
Weyl spinor of a free Majorana particle obeys the wave 
equation (see, e.g., [ZJj), 



i?)a + (cp)??a + im a a 2 r)* = 0. 



(8.3) 



Note that Eq. (|8.3p can be formally derived from Eq. (|2.7[) 
if impose the Majorana condition ip a = {ip a ) c , where the 
index "c" means the charge conjugation, on the four com- 
ponent spinor ip a . The Majorana condition means that 
this spinor is represented as ipj = (i^Vai Va)- 

The general solution of Eq. (I8.3[) can be presented 
as Fri 



x £ [a^(p)^)(p)e- i 
C=±i 



E„t 



■a^ m (-p)v^(- P y E ^], (8.4) 



where E a — <J m 2 + |p| 2 . As in Sec. [TT| the coefficient 
p) is time independent and its value is determined 



by the initial condition (|8.2 

Th 
form 



The basis spinors uffl and vi^ in Eq. (|8.4I) have the 



Wo(p)=A w_, u~(p) = A - 



E a 



E a + \p : 

where w± are helicity amplitudes given by [7 



10+ , 



(8.5) 



u+ 



,-»*/2 cos (tf/ 2 ) \ 
e^/ 2 sin(i?/2) i 

e #/2 cos (,9/2) 



the angles <fi and i? giving the direction of the momentum 
of the particle, p = |p| x (sin ■& cos <p, sin $ sin (f>, cos $) . 
The normalization factor A a in Eq. ()8.5|) can be chosen 
as 



K = 1 — 



(£ a + |p|) 2 ' 



(8.7) 



Let us mention the following properties of the helicity 
amplitudes tu±: 



(crp)w ± = ±|p|w±, icr 2 w± = +w T , u?±(-p) = iui T (p) 
(w + ® tu^) — (w_ ® 10+) = icr 2 , ^w+ <g> + flO_ ® toLj 



I 

which can be immediately obtained from Eq. (18.61) and eigenstate: 
which are useful in deriving the results given below. 

The time- independent coefficients a„(p) in Eq. (|8.4[) 
have the following form (75|: 



<(P) 



a a (P) 



'(2tt) 3 / 2 



r?i 0)t (- P K + (p) 



+ 7F 1 Tb-a +t (-p)^ 0) (p) 
Ea + |P| 



1 



(2tt) 3 / 2 
im 



u^(p)»?£ 0) (p) 



^>t(-p) u -(- p ) 



(8.9) 



where (p) is the Fourier transform of the initial wave 
function r] a , 



Va 0) (p) = J d 3 pe- ipr r?( 0) (r). 



Using Eqs. 



we then obtain the following ex- 



Va(r,t) = 



(2nY 



e lpr A 



„-i-E„t 



E„ 



iE„t 



X ® 



\E„t 



\E„t 



E a + \p 
X (w+®wty Ua^iP) 



2-^-sMEat)a 2 ^>{^) 

Ea + \P\ 



(8.10) 



From Eqs. 



and (|8.10l) it follows that a mass 



pression for the wave function for the neutrino mass eigenstate particle initially in the left polarized state 
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Tji ° (r) ~ ui_(k)e lkr is described at later times by 



{ «■*- 


m a 




E a + \k\\ 



iE a t I „ik 



e lkr w_(k) 



2i- 



sin(£; a t)e- ikr w + (k) \. (8.11) 



E a + |k| 



Let us notice that the second term in Eq. (|8.11[) describes 
an antineutrino state. Indeed the spinor u> + (k) satis- 
fies the relation, (erk)u; + (k) = |k|u;+(k), see Eq. (18. 8|) . 
Therefore it corresponds to an antiparticle (see Ref. |77|). 
This term is responsible for the neutrino-to-antineutrino 
flavor state transition v\ ■<-» (^) c . 

According to Eqs. (|2.4[) and (18. ip . the wave function 
of the left-handed neutrino of the flavor "a" is = 
cos 6r)\ - sin Orfe. From Eqs. (|B~Tj) and (|5TT|) it then 
follows that the probability of the transition vk 
vacuum is given by 



m 



■ cos(at) sin($t) 



2|k| 

x [m 2 sin(Eit) — m\ sin(i?2£)] 

,4 



Ikl 4 



(8.12) 



where a and $ are introduced in Eq. (|2.21[) . 

We can compare Eq. (|8.12l) with analogous expression 
for Dirac particles (|2.20[) . The leading terms in both 
equations coincide, whereas the next-to-leading terms 
have different signs. Therefore one can reveal the na- 
ture of neutrinos, i.e. say if neutrino mass eigenstates are 
Dirac or Majorana particles, studying the corrections to 
the transition probability formulas. 

Analogously we can calculate the transition probability 
for the process v\ — > (z>„) c using the second term in 
Eq. (EH]), 



p„ 



=I(^)T 



sin 2 (20) 
4|k| 2 



mi sin(Eit) — 7712 sin(E2t)Y 



O 



Ikl 4 



.13) 



Note that the next-to-leading term in Eq. (|8.12l) and lead- 
ing term in Eq. (|8.13|) have the same order of magnitude 
~rnl/\k\\ 

Before moving to consider Majorana neutrinos in mag- 
netic fields in Sec. IIX1 we make a general comment con- 
cerning the validity of our approach based on relativistic 
field theory involving (classical) first quantized Majorana 



fields. It has been stated [78j that the dynamics of mas- 
sive Majorana fields cannot be described within the clas- 
sical field theory approach due to the fact that the mass 
term of the Lagrangian, 77 T icr 2 /7, vanishes when 77 is rep- 
resented as a c- number function (see, e.g., Eq. (|11.27p 
below). Note that Eq. (|8.3p is a direct consequence of 
the Dirac equation if we suggest that the four-component 
wave function satisfies the Majorana condition. There- 
fore a solution of Eq. (|8.3p . i.e., wave functions and en- 
ergy levels, in principle does not depend on the existence 
of a Lagrangian resulting in this equation. The wave 
equations describing elementary particles should follow 
from the quantum field theory principles. However quite 
often these quantum equations allow classical solutions 
(see Ref. [lH for many interesting exam ples ) . We demon- 
strate in Sees. UK] (see also Refs. (H^H [M S3|) that 
oscillations of Dirac neutrinos in vacuum and various ex- 
ternal fields can be described in the framework of the 
classical field theory. The main result of this section was 
to show that the quantum Eq. (|8.3p for massive Majorana 
particles can be solved [see Eq. (|8.1ip ] in the framework 
of the classical field theory as well. 



IX. MAJORANA NEUTRINOS IN MATTER 
AND TRANSVERSAL MAGNETIC FIELD 

In this sections we apply the formalism developed in 
Sec. I Villi to the more general case of Majorana neutrinos 
propagating in background matter and interacting with 
an external magnetic field. We show that the initial con- 
dition problem of the mass eigenstates can be expressed 
in terms of the Schrodinger equation with an effective 
Hamiltonian which coincides with previously proposed 
one H3- 

For describing the evolution of two Majorana mass 
eigenstates in matter under the influence of an external 
magnetic field, the wave equation (|8.3p is to be modified 
to the following form: 

\T] a + [ap - —j T] a + im a a 2 r] a 

- -jVb - in(crB)<j 2 eabVb = °j a ¥= b > (9- 1 ) 

where e a b — i( cr 2)a6, and g a and g were defined in 
Eq. (|3.5p . Note that Eq. (|9.ip can be formally derived 
from Eq. (|5.3p if one neglects vector current interactions, 
i.e., replace (1— 7 5 )/2 with — 7 5 /2, and takes into account 
the fact that the magnetic moment matrix of Majorana 
neutrinos is antisymmetric (see, e.g., Ref. 80]). We will 
apply the same initial condition (|8.2p as in the vacuum 
case. It should be mentioned that the evolution of Ma- 
jorana neutrinos in matter and in a magnetic field has 
been previously discussed in Ref. (8lj |. 

The general solution of Eq. (|9.ip can be expressed in 
the following form: 
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'.(r,t) - / 7^72e ipr E [^ ) (p,t)^ ) (p)exp(-i40t) + ap(-p,t) U P(-p)exp(i£Pt)], (9.2) 



C=±i 



where the energy levels are given in Eq. (|3.7|) (see Ref. (HI). The basis spinors in Eq. (|9.2[) can be chosen as 



u a (P) = A a w-, u+(p) = -A+ 



K + + (|p| - ffa/2) 



^(p) = A+w_, u a (p) = A a 



E- + (\p\+g a /2) 



'+• 



where the normalization factors £ = ±, are given by 

(A«r 2 = i 



[Ea + (| P | - Cffa/2)] 



2 ' 



(9.3) 



(9.4) 



Let us consider the propagation of Majorana neutrinos in the transversal magnetic field. Using a similar technique 
as in the Dirac case in Sec. [V]and assuming k 3> m a , we end up with the following ordinary differential equations for 
the coefficients 



i^f - = H'V, H 
dt 



f 


ge iul - t /2 


fiBe in - f 


ge~ iul - t /2 





-fiBe-^+t 







-ge iul + t /2 







_ ge -i u+ t/ 2 o 



(9.5) 



where \fr T = (a x , a 2 , af, aj) and 



w± = - E 1 * « 2$ =F 



5i - 32 



n T = Ef - Et « 2$ ± 



.91 + .92 



which should be compared with Eq. (15.71) 
By making the matrix transformation 



*' =U% U = diag | e i(*+si/2)t ) e -i(*-»/2) tj e i(»- fll /2) tj e -i(* +fla /2)t j ; 
we can recast Eq. (|9.5[) into the form 

/$+5i/2 3 /2 M S \ 

-$ + .g 2 /2 
-/iB $-3i/2 -g/2 

V MS - 5 /2 -$-g 2 /2j 



(9.6) 



(9.7) 



dt 



(9.8) 



Let us note that the analogous effective Hamiltonian 
has been used in describing the spin flavor oscillations 
of Majorana neutrinos within the quantum mechani- 
cal approach (see, e.g., Ref. (20[) if we use the basis 

Note that the consistent derivation of the master 
Eq. (|9.f p should be done in the framework of the quan- 
tum field theory (see, e.g., Ref. |78j). supposing that the 
spinors r\ a are expressed via anticommuting operators. 
This quantum field theory treatment is important to ex- 
plain the asymmetry of the magnetic moment matrix. 
However, it is possible to see that the main Eq. (|9.f I) can 
also be reduced to the standard Schrodinger evolution 
Eq. (|9.8p for neutrino spin flavor oscillations if we suppose 
that the wave functions rj a are c-number objects. That is 
why one can again conclude that classical and quantum 



field theory methods for studying Majorana neutrinos' 
propagation in external fields are equivalent. 



X. SPIN FLAVOR OSCILLATIONS OF 
MAJORANA NEUTRINOS IN THE EXPANDING 
ENVELOPE OF A SUPERNOVA 

In this section we discuss the application of the formal- 
ism developed in Sec. HXI to the description of spin flavor 
oscillations of Majorana neutrinos in an expanding en- 
velope of a supernova and compare the results with the 
case of Dirac neutrinos studied in Sec. I VII 

The dynamics of the system of two Majorana neutri- 
nos in matter under the influence of an external magnetic 
field is governed by the Schrodinger equation (19.81) . How- 
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ever, as in Sec. I VII the explicit analytical solution of this 
equation is quite cumbersome. That is why we again 
consider the situation when which 
results in g\ = — gi (see Eqs. p. 141) and (|3.15p as well 
as Sec. IVip . In this case the eigenvalues of the Hamilto- 
nian ()9.8[) A = ±£± are given by 



1 



£± =- yJV 2 + 4(^B) 2 + 4$ 2 ± WR : 
R= y / (<f>cos28) 2 + (fiB) 2 , 



(10.1) 



where V = Gpn/V2 as in Sec IVII The time evolution of 
the wave function is described by the formula, 



*(*)= E [(UcSUfyexpi-iSct) 



C=±i 
+ 



(v ( ® V^j exp (i£ c t) 



(10.2) 



where Uq and are the eigenvectors of the Hamilto- 
nian (|9.8|) . given as 



\-*<J 



f-yc\ 

x c 

V i J 



(10.3) 



where 



x c ^ sin 20, 



Vi 



£ ( + <I>-V cos 26»/2 

Usin26> 
' 2(£ c + $ + V cos 26/2) 
V , 



2(pB) 2 (£ < + <I>) 



E c = [2£ c (£ c + $) - U 2 /2 + cos 26>] cos 29 

+ CW(£ c + $- Vcos26>/2). (10.4) 
The normalization coefficient iV^ in Eq. (|10.3p is given 



by N c = s Jl + x 2 c +y 2 c + z 2 . 

Proceeding along the same lines as in Sec. IVII we ob- 
tain from Eqs. ([87T]) and (|T072|) - (|T0T4|) the probability of 
the process v\ — > as, 

P v L^ uR (t) =[C+ cos{£+t) + C- cos(£_<)] 2 



where 



[S+ sm(£ + t) + SL sin(£_t)] 2 , (10.5) 



—2 [sin26»(a; c + y c z c ) + cos2% c - x c z c )}, 



1 



s < =j^(yc + x c z c)- 

Consistently with Eq. 
wave function as 



(10.6) 

we have taken the initial 



With help of Eqs. (TICU)) and ((106)) it is easy to check 
that C+ + C- = guaranteeing P(0) = 0. 

Note that formally Eq. (|10.5j) corresponds to the tran- 
sitions v\ — > u^. However, virtually it describes os- 
cillations between active neutrinos v\ «-» (f^) c since 
= f° r Majorana particles. 

As in the case of Eq. (16.7p . Eq. (jl0.5[) can be treated 
analytically for relatively small values of the effective 
potential V. The ensuing envelope functions P u d — 
Pq ± yj P 2 + P 2 depend on the coefficients and 
in the same way as in Eq. (|6.10p . The transition prob- 
abilities at various values of the matter density and the 
magnetic field are presented in Fig. [5] 

Despite the formal similarity between Dirac and Ma- 
jorana transition probabilities [see Eqs. (16.71) and (|10.5j) ] 
the actual dynamics is quite different in these two cases, 
as one can see by comparing Figs. [H[d)-|l|d) and Fig. [SJ 
panels (b), (d) and (f). In particular, in the Majorana 

case pi max - ) = 4|C+C_| for arbitrary B, to be compared 
with Eq. (I6.12p . while the function p^ max ^ has the same 
form as in the Dirac case given in Eq. f|6 . 1 3[) . In contrast 
with the Dirac case, the averaged transition probability 
does not achieve its maximal value at some moderate 
magnetic field P max value, but both _p™ ax and P™ ax are 
monotonically increasing functions of the strength of the 
magnetic field with 1 and 1/2 as their asymptotic values, 
respectively. We can understand this behavior when we 
recall that, at [iB ^> max($, U), the effective Hamilto- 
nian Eq. (|9.8p becomes 



Hoo = i/.iBj , i7 









0-10 
-10 
/ 



(10. 



(sin 6, cos 0, 0,0). 



(10.7) 



The Schrodinger equation with the effective Hamilto- 
nian (110.8)) has the formal solution 

= cxp(-ii/ oo i)*(0) = [cos(fiBt)+~f 2 sm(fiBt)}^(0). 

(10.9) 

Using Eqs. (JHTTJ) , (111377)) and ftUM we then immedi- 
ately arrive to the following expression for the transition 
probability, P(t) = \v^\ 2 = sin 2 (fiBt), which explains 
the behavior of the function p/ l max - ) at strong magnetic 
fields. Note that the analogous result was also obtained 
in Ref. [59|. 

Finally, it is worth of noticing that in contrast to the 
Dirac case, the behavior of the transition probability in 
the Majorana case is qualitatively similar for different 
matter densities and different magnetic fields [see Fig. [SJ 
panels (a), (c), and (e)]. 

The problem of Majorana neutrinos' spin flavor oscil- 
lations was studied in Refs. [82|, I83I with help of numer- 
ical codes. For example, in Ref. |82j the evolution equa- 
tion for three neutrino flavors propagating inside a pre- 
supernova star with zero metallicity, e.g., corresponding 
to W02Z model (58[, was solved for the realistic mat- 
ter and magnetic field profiles. Although our analytical 
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x, 10 3 km [IB, 10 11 eV 



(e) (f) 

FIG. 5. (color online) (a), (c) and (e) The transition probability versus the distance passed by a neutrino beam. The neutrino 
parameters have the same values as in Figs. OH E v — 10 MeV, 8m 2 — 8 x 10~ 5 eV 2 and 8 — 0.6. The magnetic energy is equal 
to \iB — 10 -12 eV. The black line is the function P(x), the blue and green lines are the envelope functions P u ,d(x) and the red 
line is the averaged transition probability Po(x). (b), (d) and (f) The dependence of the maximal values of the functions P(x) 
and Po(x), blue and red lines respectively, on the magnetic energy for the given density. Panels (a) and (b) correspond to the 
matter density p — lOg/cc, (c) and (d) to p — 50g/cc and (e) and (f) to p = lOOg/cc. This figure is taken from Ref. [3S|. 



transition probability formula (|10.5I) is valid only for the 
constant matter density and magnetic field strength, it 
is interesting to compare our results with the numeri- 
cal simulations of Ref. 821. In those calculations the 



authors used magnetic fields B ~ 10 10 G and magnetic 
moments ~ 10~ 12 /xb that give us the magnetic energy 
liB r~j 10~ n eV. This value is the maximal magnetic 
energy used in our work. 
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It was found in Ref. [82| that spin flavor conversion 
is practically adiabatic for low-energy neutrinos corre- 
sponding to E v ru 5 MeV inside the region where Y e sa 0.5 
and the averaged transition probability for the channel 
— > v e is close 0.5. This big transition probability is 
due to the RSF-H and RSF-L resonances at the distance 
i=3 O.Oli?©. Even so that we study two neutrino oscil- 
lations scheme, we obtained the analogous behavior of 
p(max) j. gee ^ However, in our case this big tran- 
sition probability is due to the presence of the strong 
magnetic field (see Refs. [3 HI])- We cannot compare 
our transition probability formula (|X0.5|) with the results 
of Ref. [HJ for higher energies, E v > 25 MeV, since spin 
flavor oscillations become strongly nonadiabatic for these 
kinds of energies and one has to take into account the co- 
ordinate dependence of the matter density which should 
decrease with radius as 1 /r 3 (84j • 



XI. EVOLUTION OF NEUTRINOS EMITTED 
BY CLASSICAL SOURCES 

In this section we present an alternative approach to 
the description of neutrino flavor oscillations which in- 
volves the studies of the evolution of mixed massive neu- 
trinos in vacuum under the influence of external classical 
fields [38j |. These fields are supposed to be localized in 
space and treated as "sources" which emit flavor neutri- 
nos. 

As in Sees. ITlWl we consider a number of flavor neu- 
trinos which, in principle, can be arbitrary: A = 
e, fj,,T, . . . . Note that we will formulate the dynamics 
of the system in terms of the left v\ and right handed 
chiral components of the neutrino spinors. 

The external sources are denoted as J A . It should be 
noted that, if one describes the emission of mixed mas- 
sive neutrinos in any process within the quantum field 
theory, there is a nonzero probability of emission of neu- 
trinos of different flavors (85( . Although a source of one 
particular flavor can prevail, one should always consider 
the situation when the sources of all flavors are present. 

The Lorentz invariant Lagrangian describing the evo- 
lution of this system has the form, 

c = E $?f d A + ^vax) 

A 

AA' 

+ ™AV KfC^+h.C.) 

+ E(^VA^ + h.c), (H.l) 

A 

where (m^ x ,), (wi^A')' ( m AA') are the Dirac as well as 
(left and right) Majorana mass matrices and C = ij 2 -f° 
is the charge conjugation matrix. These mass matrices 
should satisfy certain requirements which are discussed 



in Ref. [44j . The Lagrangian (ITTTTj) is CPT- invariant. In 
this section we adopt the notations for Dirac matrices as 
in Ref. . The analogous mass terms are generated in 
theoretical models based on the type-II seesaw mecha- 
nism (see, e.g., Ref. |45|). 

The external fields J A = J A (r, t) can be arbitrary func- 
tions. If we suppose that neutrinos interact with external 
sources in frames of the electroweak model, the spinor £\ 
in Eq. (|ll.ll) is the charged SU(2) isodoublet partner of 
V\. For example, we can study the neutrino emission in 
a process like the inverse /3-decay: p + e~ — > v e + n. In 
this case the external fields in Eq. are (see, e.g., 

Ref. |3) 

J£ = - V2G F $ n ^(l - a 7 5 )^, 

J£=0, J£=0, (11.2) 

where Hf p and & n are the wave functions of a proton and 
a neutron, a « 1.25. In Eq. (|11.2[) we assume that a 
source consists of electrons, protons and neutrons. 

To proceed in the analysis of the dynamics of the sys- 
tem (jll.ip we should diagonalize the mass term of the 
Lagrangian. The diagonalized Lagrangian is expressed in 
terms of the Majorana fields with different masses [IH. In 
general case the number of these Majorana fields is dou- 
ble than the number of flavors in the initial Lagrangian. 
In the following we discuss the case when only Dirac mass 
term is presented. Then we study the general situation. 

A. Dirac mass term 

In this section we suppose that only Dirac mass matrix 
is present in Eq. (jll.ip . Analogously to the discussion 
in Sec. [TT] we introduce the new set of spinor fields %j) a . 
In our situation the mixing matrix is a square unitary 
matrix (t/ Aa ), 

y\ = ^Uxai> a . (n.3) 

a 

By definition the mass eigenstates ip a are Dirac particles. 

When we transform the Lagrangian using 
Eq. (|11.3|) . it is rewritten in the following way: 

C = M^d,, - m Q )V> Q + E$»& + h.c), (11.4) 

a a 

where £ a is the external source for the fermion ifi a , 

e^E^A^A/r- (n.5) 

A 

Using Eqs. (|11.4I) and (|11.5p we receive the inhomoge- 
neous Dirac equation for the fermion tp a , 

(i7"0 M - m„W„ = (11.6) 

As in Sec. [TT] the masses m a are the eigenvalues of the 
matrix (m AA ,). 
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The solution to Eq. (|11.5I) for the arbitrary spinor £ a 
is expressed with help of the retarded Green function for 
a spinor field (see Ref. |88j|). 

^.(r,t) = J d?v'At'ST{v-v',t-t')U^X)- (H.7) 
The explicit form of S^ ct (r, t) can be also found in 



Ref. [gal, 

^ et (r,i) = (i 7 ^^+m a )^ et (r,t). (11. 



In Eq. (|11.8JI _D' ct (r, i) is the retarded Green function for 
a scalar field [8811 , 



d 4 p 



(27r) 4 m.2 - p 2 + iep° 

LewU^-Otf^J^maa)}, (11.9) 
In I 2s J 



where 8(t) is the Heaviside step function, J\{z) is the 
first order Bessel function, s = \Jt 2 — r 2 and r = |r|. 

To proceed in further calculations it is necessary to 
define the behavior of the external sources. We assume 
that the sources are localized in space and emit harmonic 
radiation, 

t x (r,t)j£(T,t) = 9(t)l x J^e- im S 3 (r), (11.10) 

where l\ = is the time independent component of the 

spinor i\, J\° is the amplitude of the function Jf(r, t), 
and E is the frequency of the source. Using Eqs. (|11.5[) 
and (jll.lOp we obtain for £ n , 



e^E^M ^ (ii.il) 



With help of Eqs. (fTTTH|) - (j 1 1 . 1 1 j> we can rewrite Eq. (fTTTT]) 
in the following way: 

/ d 
^ a (r,t) = ^—+m a 

~ iEt J*dTD™\r,T)e iE A&\ (11.12) 

One can see that two major terms appear in Eq. (|11.12[) 
[see also Eq. (| 1 1 .9(1 ] while integrating over r 



x < e 



drD™ t (r,r)e iBT =I 1 +I 2 . 



They are 



2tt 



6{t - r) - Er 



(11.13) 



(11.14) 



| Ji(m a s) ciEr 



Ji{m a x) ciE ^pi+^i ^ u lg s 
\/r 2 + x 2 



where x m = \Jt 2 — r 2 and s = \Jt 2 — r 2 . It is interesting 
to notice that both X\ and Z 2 in Eqs. (|11.14[) and (|1 1 . 15[) 
are equal to zero if t < r. It means that the initial 
perturbation from a source placed at the point r = 
reaches a detector placed at the point r only after the 
time t > r. 

Despite the integral X\ in Eq. (|11.14[) is expressed in 
terms of the elementary functions, the integral Z 2 in 
Eq. (|11.15[) cannot be computed analytically for arbi- 
trary r and t. However some reasonable assumptions can 
be made to simplify the considered expression. Suppose 
that an observer is at the fixed distance from a source. As 
we have already mentioned one detects a signal starting 
from t > r. It is obvious that a nonstationary rapidly 
oscillating signal is detected when the wave front just 
arrives to a detector, i.e. when t > r. The situation 
is analogous to waves propagating on the water surface. 
Therefore, if we suppose that one starts observing parti- 
cles when the non-stationary signal attenuates, i.e. t 3> r 
or x m — > oo, we can avoid relaxation phenomena. In this 
case the integral I 2 can be computed analytically, 



1 

Anr 



(e 



iEr 



), 



(11.16) 



where p a = y/ E 2 — m 2 is the analog of the particle mo- 
mentum. 

To obtain Eq. (|11.16l) we use the known values of the 
integrals, 



da;- 



J ^ mx) sin(EV^T^) 



yjr 2 + 

7T 

= 2 Jv/2 



x J_ 



v/2 



E + VE 2 - 



dx- 



J v (mx) cos ^^/ r 2 _|_ x 2^ 



yjr 2 + 

7T 



- [E-y/E?-m 2 



x N. 



-i//2 



t -[e + VW~ 



(11.17) 



and the fact that the Bessel and Neumann functions of 
±1/2 order 



(11.18) 



Ju 2 (z) =N_ 1/2 (z) = d — sinz, 



J-l/2(z) =\ COS 



are expressed in terms of the elementary functions. The 
approximations made in derivation of Eq. (|1 1 . 16[) are 
analysed in Appendix [D] 

Using Eqs. (|11.12|) - (lll.l8j) we obtain the field distribu- 
tion of the fermion %p a in the following form: 



Mr,t) = e- iEt+ip " r O a 



AO) 

So 

Airr ' 



(11.19) 
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where O a — j°E — (-fn)p a + m a and n is the unit vector 
towards a detector. It should be noted that in deriving 
Eq. (|11.19j) we differentiate only exponential rather than 
the factor 1/r because the derivative of 1/r is propor- 
tional to 1/r 2 . Such a term is negligible at large distances 
from a source. We also remind that Eq. (|11.19|) is valid 
for t 3> r. 

Let us turn to the description of the evolution of the 
fields v\. Using Eqs. (|11.3[> and (|11.19j) . we obtain the 
corresponding wave function, 



fA(r,t) 



-\Et 



a A' 



t(0)m 

Aixr 



(11.20) 



As in Sees. InWl we study the evolution of only two 
fermions. The mixing matrix is given in Eq. (|2.4[) . We 
also choose the amplitudes of the sources in the following 



way: 



0, J, 



This choice 



= J»^0 and If) 
of the sources corresponds to the emission of neutrinos of 
the flavor "/3" and detection of the neutrinos belonging 
to the flavor "a" at the distance r from a source. Using 
Eqs. (|2.4I) and (111.20[) we obtain the wave functions of 
each of the particles y a ,pi 



J 11 



v a (r, t) = sin 9 cos Be lEt — — 
Airr 

x (e^Oi - e ip2r 2 )7^, 



V/3(r,t) 



,-lEi 



Anr 



(sin 2 9e ipir 1 + cos 2 0e ip2r O 2 ) 7 ^. (11.21) 



In the following we discuss the high frequency approxi- 
mation, E 3> mi 2i which corresponds to the emission of 
ultrarelativistic neutrinos. 

The probability to detect a neutrino of the flavor "A" 
can be calculated as P\(r,t) = \v\(r,t)\ 2 . Finally using 
Eq. (|11.21j) we get for the probabilities, 



/tM J" 

P p {r)=-2E 2 j—^{T, vl 



(47T7-) 



x <^ 1 -skr(2(9)sin 2 



8m 
~A~E 



O 



E 2 



(11.22) 



where (T^ v ) — l*^J l (aaa)') u l, It is possible to calculate 
the components of the tensor (T„ y ). They depend on the 
properties of the fermion £p = £, 

(Too) = -H, (Toi) = (Tfl,)* =rn- i[n x C]i, 
(Tij) = iy(vn) - (vifij + VjUi) + \e ijk n k , (11.23) 



where v = (a) is the velocity of fermion £ and C = (S) 
is its spin. 

Let us discuss the simplified case when the spatial com- 
ponents of the four-vector J' 1 arc equal to zero: J = 0. 
It corresponds to the neutrino emission by a nonmoving 
and unpolarized source. Using Eqs. (|11.22j) and (|11.23p 
we can find the transition and survival probabilities in 
the following way: 



Pp->a(r) ~sm 2 (26>)sin 2 



AE 



Pp^fi{r) -1- sin 2 (26») sin 2 



Sm 2 
— — ) 
AE 



(11.24) 
(11.25) 



where we drop the common factor 2i? 2 (vn)| J°| 2 /(47rr) 2 
to have the probabilities normalized to unity. It should 
be noted that Eqs. (I11.24[) and (|11.25l) are the same as the 
common formulae for the description of neutrino flavor 
oscillations in vacuum HI, [43 . 



B. General mass term 

When both Dirac and Majorana mass matrices are 
present in Eq. (|11.1[) . left and right handed chiral com- 
ponents should be transformed independently (see also 
Ref. 0), 



E^c 



E 



(11.26) 



with help of the matrices (U\ a ) and (V\ a ) to diagonal- 
ize the general mass term. Note that these matrices are 
rectangular and nonhermitian. The modern parameteri- 
zation for these matrices is given in Ref. [89]. 

However one cannot say whether particles f^' R are 
Majorana or Dirac although they correspond to definite 
mass eigenstates. That is why we introduce a new field 
%jj a = ^ + ( v I / a) c which is Majorana by definition. 

The Lagrangian for the particles ip a takes the form, 



C = 



(V^^a - m,^ 



(11.27) 



where the external sources £ a = £j have the same form 

as in Eq. (|11.5j) . 

As we mentioned in Sec. I Villi Majorana spinors are 
equivalent to two component Weil spinors [75[. Hence 
we can rewrite the spinors f/'a' R an d £a as 



where 



^E^Wa - 



(11.29) 
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To obtain Eq. (jll.29|) we suppose that, as in Sec. IXI A[ 
the vectors j£ have only time component J°. We also 

assume that (i\) T = (0,xa)- To derive Eq. (|11.29|) it 
is crucial that only left handed currents interactions are 
presented in Eq. ([ll.ip . 

It is useful to rewrite the Lagrangian (|11.27[) in terms 
of the two component spinors r] a and (fi a [74j | , 



a 

+ E ( 2 m<i7? « ' 2r? « + ^ a + h - c ' 



(11.30) 



Using Eq. (|11.30[) we can receive the wave equation for 
two-component spinors, 



— - erV ) rj a + m a a 2 r\ a = \(p a 



(11.31) 



The solutions of Eq. (|11.31|) have the form (see, e.g., 
Ref. Ml 

Va(r,t) = / - d 3 r / dt'^ et (r-r',t-t / )0a(r , ,t'): (U-32) 
r,* a (r,t) = ydVdtX ot (r-r',<-i')0a(r',i'), (H-33) 

where the retarded Green functions are expressed as (see 
also Ref. Q), 



S' et (r,t) =ic^Df (r,t), 



(11.34) 



Here <9 M = (<9 t ,V), a* = (a°,tr), a is the 2 x 2 unit 
matrix, <x are Pauli matrices, and D„ et (r, f) is given in 
Eq. (flOj) . One can check that Eq. (|11.32p and (|11.33j> . 
along with the definition of the retarded Green func- 
tions (|11.34l) . represent the solutions to Eq. (|1 1 .31[) by 
means of direct substituting. 

Let us assume that the sources 4> a depend on time and 
spatial coordinates as in Sec. IXI Al 



^„(r,t) =e(t)^e- iEt 5 3 (r), 

^ 0) =E^W1 0) , 



(11.35) 



where X\ is the time independent component of \\ ■ in 
deriving of Eq. (|11.35j) from Eq. (|11.29|) we introduce 



the new quantities J\ = , where J\''" is the time 
independent part of J°, to simplify the notations. 

On the basis of Eqs. (|11.32p - (|11.35p and using the tech- 
nique developed in Sec. IXI Al we get the particles wave 
functions, 

,(o) 

Va (r,t) =i^d^e- iEt+ ^ r }^, (11.36) 

47rr 

,(o) 

v * a (r,t)=im af j 2 e- iEt +^ r ^-. (11.37) 



r(0)0 



To derive Eqs. (I11.36[) and (|11.37[) we suppose that t > r. 
The derivatives in Eq. (|11.36[) are applied on the expo- 
nent only because of the same reasons as in Sec. IXI Al 

Using Eqs. (|Ti~26l) . (jTi~36l) and (|Ti~37l) as well as the 
following identity: 

(^) C =K) R =E^X, 

a 

("x) C =cfi) T , (H.38) 
we get the wave functions of v\ and {yY) c as 

aA' 

(11.39) 




W)°('.*)=-^e- UM E^ r ^ 

aA 7 




(11.40) 



where x^J? = (1/2) [1 — (cm)]xy' ■ I n deriving of 
Eq. (|11.39j) we suppose that 



,(0) 



[E-p a (trn)]*E[l-(<m)], 



(11.41) 



that is valid for relativistic neutrinos. 

It is possible to construct two four component Majo- 
rana spinors from two-component Weil spinors, 



(i) 




(2) 




(11.42) 



where r) a and 77* are defined in Eqs. (|11.36|) and (|11.3T|) . 
One can see that these spinors satisfy the Majorana con- 



dition 



1.2) 



ip^' 2 \ We use the spinor ip^' to re- 
ceive Eq. (|11.39P and ip^ for Eq. (|11.40|) . In this case 
we get that v\ and {y\Y are obtained as a result of the 
evolution of particles emitted from the same source. 

It should be mentioned that both v\ and \V\\ in 
Eqs. (|11.39[) and (|11.40[) propagate forward in time. To 
explain this fact let us discuss the complex conjugated 
equation (|11.32p . Performing the same computations one 
arrives to the analog of Eq. (|11.40p which would depend 
on time as e lEt . However this wave function describes 
a particle emitted by the source different from that dis- 
cussed here. Indeed, if we studied the complex conju- 
gated Eq. (|11.32[) , the integrand there would b e S™**^*. 
It would mean that the source in Eq. (I11.35P would be 



proportional to 



1.(0)* 



\Et 



that, in its turn, would sig- 



nify that {yY) c would be emitted in a process involving a 
lepton which is a charge conjugated counterpart to that 
discussed here. 
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Let us illustrate this problem on a more physical ex- 
ample. Suppose that we study a neutrino emission in a 
process like the inverse /3-decay, 

p + q -^n + vx^n + ^Uxa^, (11-43) 

a 

where p, n and £ x stand for a proton, a neutron and for 
a negatively charged lepton. The complex conjugated 
Eq. (|11.32[) would correspond to a process, 

n+e+^ P + i>x^p + J2 u L W) c 

a 

a 

where % and l\ denote an antineutrino and a positively 
charged lepton, which is the charge conjugated counter- 
part to £ x . In Eq. (|11.44j) we use the facts that only left 
handed interactions exist in nature and the fields ip a de- 
scribe Majorana particles. As one can see, Eqs. (|11.43[) 
and (|11.44[) represent two different processes. Hence, if 
we used (rj a )* to obtain (r,t), i.e. after a beam 

of neutrinos passes some distance r, it would correspond 
to the initial reaction (|1 1.441) rather than (|11.43|) . Note 
that the same result also follows from Eq. (|11.42p if we 
replace ipcP «-» tp^ there. 

Let us suppose for simplicity that the momentum of 
the fermion l\ is parallel to the neutrino momentum. It 
takes place if a relativistic incoming lepton is studied. 
We also assume that this fermion is in a state with the 
definite helicity, 

i[l-(<xn)]xi 0) =xl 0) . (11.45) 

This expression is again natural for the relativistic 
fermion l\. One can notice that the Lagrangian (|11.1[) 
is written in terms of the left handed chiral projections 
of l\. Therefore, if we study a relativistic lepton, it will 
have its spin directed oppositely to the particle momen- 
tum as one can see from Eq. (|11.45[) . 

Using Eqs. (fTl~39l) . (fl~L40| and (fTPS) as well as 
the orthonormality of the two-component spinors xl i 
(x^Xy^ = <5aa', we get the probabilities to detect v\ 
and [y\Y as 

X U Xa Ul y UlUy b \Jy\ 2 , (11.46) 
P, ,c(r\ ~ V maTOfc r i (Pa-p t .)r 

K) V > ^ (2E) 2 

ab,\' y ' 

*Ui a Ul x ,Uj x U Wb \Jy\ 2 . (11.47) 

In Eqs. (|11.46l) and (|11.47|) we drop the factor 
(2£)7(47rr) 2 . 



We can see that Eqs. (|11.46|) and (|11.47[) contain the 
oscillating exponent. This our result reproduces the 
usual formulae for neutrino oscillations in vacuum (44j. 
It should be also noticed that the expressions for the 
probabilities depend on r rather than on t in contrast to 
our previous works (HHIt], HH, H(| ■ Note that the prob- 
lem whether neutrino oscillations happen in space or in 
time was also discussed in Ref. [90]. The similar coor- 
dinate dependence of the probabilities was obtained in 
Refs. [26l |27| where the problem of neutrino oscillations 
in vacuum was studied. 

It was mentioned in Ref. [26| that oscillations between 
active and sterile neutrinos are possible in case of the 
nonunitary matrix (U\ a ), i.e. the presence of only Ma- 
jorana mass terms is not sufficient for the existence of 
this kind of transitions. The situation is analogous to 
that considered in the pioneering work [911 ] where oscilla- 
tions between neutrinos and antineutrinos were studied. 
In Eq. (|11.47[) we obtain that the probability to detect 
{yY) i s suppressed by the factor m„,rrih I ' E 2 . It is in 
agreement with the results of Refs. [2g, [92| as well as 
with Eq. (|5T3l) . 



XII. QUANTUM FIELD THEORY 
DESCRIPTION OF NEUTRINO OSCILLATIONS 
IN BACKGROUND MATTER 

In the present section we study neutrino oscillations in 
background matter (see Sec. IIIII) in frames of the quan- 
tum field theory approach. In particular we will be in- 
terested to examine the influence of background matter 
on the transitions between neutrinos and antineutrinos 
discussed in Sees. I Villi and IXJBl 

Despite we could receive satisfactory results for the 
description of neutrino-to-antineutrino transitions in vac- 
uum in frames of the relativistic quantum mechanics ap- 
proach [see Eqs. (|8.13|) and (jll.471) ] the accurate treat- 
ment of this process has to be done within the quantum 
field theory because of the following reason. For the case 
of Majorana neutrinos, particles are identical to their an- 
tiparticles and only in frames of the quantum field the- 
ory one has the most accurate description of antiparticle 
states. Moreover, as we have seen in Sec. IVIII1 "antineu- 
trino" states appear along with the small factor m a / E u . 
Typically various approaches to the description of neu- 
trino oscillations give contradictory results at this order 
of accuracy [94[. That is why we will use the quantum 
field theory treatment to capture such a tiny effect. 

As in Sec. lXIl we will study the system of flavor neutri- 
nos v\, A = e, /i, r, . . . , propagating in dense background 
matter between two spatial points, xi and X2. We sug- 
gest that emission and absorption of neutrinos is due to 
the interaction with leptons l x and heavy nucleons N. 
i.e these interactions are localized in two spatial regions: 
a "source" and a "detector". On the contrary, the neu- 
trino interaction with background matter is uniformly 
distributed along the total neutrino propagation distance 
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L = x 2 — xi. 

In general case the dynamics of the system of mixed 
massive neutrinos should be formulated in frames of the 
left v\ and right handed projections of flavor neutri- 
nos [see also Eq. 



CvY, 



X 

- Y, ( m °A'^v + m \x' iyx 

XX' 

+ K) T C^+h.c.) 

- E ^A'ttx> - W ^ + h-c) , (12.1) 



whcr 



AA' 



and 



L,R 
l AA' 



are Dirac and Majorana mass ma- 



trices defined in Sec. DO] and 



A 



A- 



(12.2) 



is the neutrino-lepton current. The effective potential 
of the neutrino interaction with background matter /(*., 
[see Eq. p.ip ] is supposed to be coordinate independent, 
whereas the nuclear current J M is localized in space. We 
will take into account the interaction with background 
matter exactly. As in Sec. IIIII here we study the general 
case and consider the matrix (f£y) to be nondiagonal. 
Note that in Eq. (|12.1[) right handed neutrinos do not 
participate in interactions with other particles, i.e. they 
are sterile. 

As we will see below, the transitions between neutrinos 
and antineutrinos manifest in the process like (N±, N 2 ) + 
17 — s- (neutrinos) — > (N[, N£) + It- It should be noted 
that in this process massive neutrinos a ppe ar as virtual 
particles rather than show up explicitly [26l |27| . The S- 
matrix element for this kind of reaction has the form (2(| , 

S = - i ( V2G F ) 2 J d A xd\j (12.3) 
x {lt l -N' 1 N' 2 \T{f{x)J^x)f{y)J v {y)}\r -N 1 N 2 ), 
I 



where N\,2 and N[ 2 are states of the initial and final nu- 
cleons. This process is schematically illustrated in Fig|6] 
Note that the propagation of internal neutrinos is shown 
by a broad line which means that it accounts for the in- 
teraction with background matter exactly. 

Now we have to diagonalize the mass matrices m^\i 
and rn^y in Eq. (| 12 . 1[) by means of the matrix trans- 
formations (lll.26[) . with the rectangular matrices (U\ a ) 
and (V\ a ) being nonunitary and independent. The am- 
plitude of the process shown in Fig.[6]is nonzero if virtual 
neutrinos correspond to Majorana fields. That is why we 
introduce a new Majorana field ip a — ^\ + (^„) c as in 

Sec. urn 

Let us rewrite the Lagrangian (|12.1|) in terms of the 




FIG. 6. The Feynmann diagram corresponding to the lepton- 
nuclei scattering, (JVi, JVa) + 17, — > (neutrinos) — > (7V{ , A^) + 
£+, with the S-matrix element given in Eq. (112.31) . 
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ab 



where (g^ b ) is the matrix of the effective potentials of the 
neutrino interaction with matter in the mass eigenstates 
basis which is defined in Eq. (|3 . 3[) . In Eq. (|12.4I) we do 
not present the expression for the sources in the explicit 
form since they are supposed to be localized in space. 

With help of Eq. (|11.26l) the T-product of the neutrino- 
lepton currents in Eq. (|12.3I) takes the form, 



(li\T{3v(x)Mv)}\W = 



Y u « a Up b u T (-p a ){rf) (o\t{[Mx)] My)}\ohiu(p ), 

ab 



(12.5) 



where u(p a ^) are the outgoing and incoming spinors with 
four momenta p a ,/3 = (Ba,j3)Pa,i3) an d V is the normal- 
ization volume. In Eq. (|12.5I) we use the representation of 
lepton states in the form of a plane wave from Ref. (7t| . 

Since the mass of the nucleons is much bigger that the 



I 

energies of leptons, we can suggest that they are at rest, 
i.e. we replace the nucleon currents 

J M (x) = <5o M <5(x- x 2 ), J^x) = 5 p6(y - Xi), (12.6) 

in Eq. ([I2~^ . 
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The T-product of Majorana neutrino fields in 
Eq. (|f2.5p can be rewritten as (see, e.g., Ref. [74]]) 



(0|T{ [Mx)] T My)}\0) = ~CS ab (x - y), 



(12.7) 



where C is the charge conjugation matrix and S ab {x — 
y) = (0\T {i/; a (x)'ip b (y)}\0) is the usual Feynmann propa- 
gator corresponding to a Dirac particle. 

To simplify the calculations we suggest that the ma- 
trix (g^ b ) in Eq. (112.4)) is close to diagonal. Under this 
assumption the propagator in Eq. (|12.7p has the diago- 
nal form: S a b(x) — S ab S a {x). The explicit form of the 
function S a (x) is given in Appendix iDl Moreover, as in 
Sees. Hill and \V \ we will study the situation of nonmoving 
and unpolarized matter, with matrix (g^ b ) having only 
the zeroth component. 

Using the results of Appendix [D] we can represent the 
5-matrix element in Eq. (|12.3[) as 



S =2n8{E a - Ep) 
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x (1-7 5 ), 



(12.9) 



and fcf 2 = (Pa±|<?a|/2) 2 , p a = \J E 2 — m 2 is the neutrino 
"momentum" , E = E a = Ep is the energy of leptons 
which is conserved since nucleons are supposed to be at 
rest. 

Let us discuss the case of high energy leptons with 
E\ 3> m\, where m\ is the lepton mass. In this case 
the biggest contribution to the matrix element of the to- 
tal process arises from the channel in which the helicity 
of leptons changes from —1/2 to 1/2. Thus the matrix 
element in Eq. (| 12 .8[) takes the form, 

£) « ^ E Em a U* aa U; a (F a ) 



x e~ iv/2 sin 







(12.10) 



where the spherical angles tp and i? fix the direction of the 
outgoing lepton momentum with respect to the incoming 
lepton momentum, which, for simplicity, is chosen to be 
directed along the vector L. 

Using Eq. (I12.9|) we can calculate the function F a in 



Eq. (|12.10p in the high energy leptons approximation, 



(F a ) =u{p a )F a u{pp) 



(12.11) 
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The details of the averaging over the lepton states in 
Eqs. (|12.10p and (|12.11|) can be found in Ref. 0. 

Finally the total cross section for the leptons scattering 
(Ni,N 2 ) + lp -> (neutrinos) -> (N[,N£) + 2+ can be 
presented in the following form: 



a(lg -> 1+) 



GpE 2 
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ab 
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where v a a are the velocities of outgoing and incoming 
leptons. If we study the most general Lagrangian (|12.1I) . 
which contains both Dirac and Majorana mass matrices, 
Eqs. (|12.1ip and (|12.12p involve transitions between neu- 
trinos and antineutrinos studied in Refs. [26|,|92j (see also 
Sees. EEED and EJ). 

Using Eqs. (|12.1ip and (|12.12p one can conclude that 
for light neutrinos with m a <C E the matter contribu- 
tion is negligible since p a ~ E — m 2 /2E ^> |<; |. How- 
ever the situation is completely different in the case of 
heavy neutrinos which are not excluded in some theoret- 
ical models [95| . If we study the neutrino propagation in 
the dense nuclear matter with \g a \ ~ eV, the "momen- 
tum" of such heavy neutrinos can be comparable with 
the effective potential of matter interaction at the appro- 
priate choice of lepton energy. Nevertheless we can see 
in Eq. (| 12 . 1 1[) that the contribution of such heavy mass 
eigenstates to neutrino-to-antineutrino transitions will be 
strongly suppressed by interaction with background mat- 
ter. Note that such a behavior of the transition proba- 
bility in presence of background matter is in agreement 
with the results of Ref. [93| where oscillations between 
neutrinos and antineutrinos were studied in frames of the 
standard quantum mechanical approach. 

In the derivation of the main results (|12.11[) 
and (|12.12p we neglected the contribution of the non- 
diagonal elements of the matrix ($„(,)• One can expect 
that this contribution can be responsible for the enhance- 
ment of the probability of neutrino-to-antineutrino tran- 
sitions as it can for neutrino flavor oscillation in matter 
(see Sec. IIIip . If one takes into account these nondiago- 
nal elements, the neutrino propagator in Eq. (112.71) also 
acquires nondiagonal entries. It is rather difficult to ana- 
lyze this kind of "nondiagonal" propagators analytically. 
Thus, the approach for the description of neutrino oscil- 
lations based on the quantum field theory does not seem 
to be applicable for the situation of big nondiagonal ele- 
ments of the matrix (g^ b )- 
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XIII. CONCLUSION 

In conclusion we mention that in the present work we 
summarized our resent achievements in the theoretical 
description of neutrino oscillations in vacuum and exter- 
nal fields. Basically we studied three approaches: rel- 
ativistic quantum mechanics approach (Sees. iHlTVl IVIII1 
and llXI) . classical external sources method (Sec. lXIl) . and 
quantum field theory treatment of neutrino oscillations 
(Sec. EH. 

The formulation of the initial condition problem (|2.2[) 
for the system of flavor neutrinos is the main feature of 
the relativistic quantum mechanics method. Then one 
looks for wave functions of flavor neutrinos at subsequent 
moments of time. In vacuum the Cauchy problem for 
neutrino mass eigenstates (|2.3I) can be solved exactly [see 
Eq. (I2.12[) ] giving us neutrino wave functions for arbitrary 
initial condition. When we study the time evolution of 
flavor neutrinos in presence of external fields within this 
formalism, we can reduce the dynamics of the system to a 
Schrodingcr like equation [see, e.g., Eqs. (|5.10|) and (|9.8|) ]. 

Note that, using this formalism, one solves not only 
the problem of the structure of the effective Hamiltonian 
but also we clearly point out what kind of space of wave 
functions this effective Hamiltonian acts in. In frames of 
the standard quantum mechanical approach [2(1 one 
studies the evolution of a quantum mechanical neutrino 
"wave function" i/ T = {y a , v $ , . . . ), but the origin of the 
components of this wave function is disguised. Using the 
relativistic quantum mechanics method one immediately 
gets the answer to this questions. If we take the coef- 
ficients tip of the decomposition of a four component 
neutrino wave function over the basis spinors [see, e.g., 
Eqs. (|5.4[) and (I9.2[) ] and then make a certain transforma- 
tion, like in Eqs. (|5.9p and (I9.7[) . of these coefficients, we 
arrive to the quantum mechanical neutrino "wave func- 
tion" . 

Note that in our approach the coefficients are o 
numbers rather than operators as in quantum field the- 
ory. It means that this method may be called as a clas- 
sical field theory treatment of neutrino oscillations since 
we do not use second quantized neutrino wave functions. 

Besides resolving the conceptual problem of neutrino 
oscillations, the relativistic quantum mechanics method 
can be applicable for the description of neutrino oscilla- 
tions not only in vacuum (Sec. [IT]) , but also in background 
matter (Sec. IIIip . spin flavor oscillations in an external 
magnetic field (Sec. IIV|) . and in the combination of back- 
ground matter and a magnetic field (Sec. [V}. Note that 
both Dirac and Majorana neutrinos can be treated within 
this formalism (see Sees. IVlfll and HX|) . 

When we study the evolution of neutrinos in an exter- 
nal field, we use exact solutions of wave equations in this 
external field. That is why our treatment is valid for arbi- 
trary strength of external fields. Moreover, the dynam- 
ics of mixed neutrinos is usually described in the mass 
eigenstates basis. In general case, external fields, like 
interaction with background matter or with an external 



electromagnetic field, are not diagonal in this basis [see 
Eqs. (|3.3[) and (|4.4I) ]. i.e. the neutrino mass eigenstates 
are not independent in presence of an external field [see, 
e.g., Eq. (|5.3I) ]. Nevertheless in frames of the relativistic 
quantum mechanics method we can treat this kind of cou- 
pled mass eigenstates and receive results which are con- 
sistent with the standard quantum mechanical approach. 

As it was mentioned above, the relativistic quantum 
mechanics method is a field theory counterpart of the 
usual quantum mechanical description of neutrino oscil- 
lation in a sense that it allows one to describe the evolu- 
tion of mixed massive neutrinos from the first principles 
and thus it throws daylight upon some unclear issues of 
the quantum mechanical approach. However, trying to 
follow quantum mechanical approach as close as possible, 
we also adopted some of its weaknesses. For example, as 
in the quantum mechanical description, in frames of our 
method neutrino "wave functions" evolve in time rather 
than in space. It is, however, a contradiction with the 
majority of the experiments where one measures neutrino 
oscillations in space rather than in time. 

Among other disadvantages of the present method we 
mention its limited applicability to description of neutri- 
nos with small initial momentum in the case of nondiag- 
onal external fields in the mass eigenstates basis. In this 
situation the ordinary differential equations for a« (i) 

and b£\t) evolution [see, e.g., Eq. (|3.9I) ] start to entan- 
gle, making possible the transitions between "positive" 
and "negative" energy states. Thus, for example, instead 
of 4 x 4 effective Schrodinger equation for spin flavor os- 
cillations one gets 8x8 differential equation, which is 
quite difficult for the analysis. 

We should also mention the fact that to get the consis- 
tency with the quantum mechanical approach, one should 
consider rather broad initial wave packet [see Eq. (|2.18l) ] 
which is rather difficult to implement in practice. In 
Sec |TT] we also demonstrated, for the case of neutrino 
evolution in vacuum, that spatially limited initial wave 
packets do not reveal flavor oscillations since they prop- 
agate as massless particles. 

As we mentioned above, the relativistic quantum me- 
chanics approach predicts neutrino oscillations in time 
rather than in space as they are observed in experiments. 
To resolve this contradiction in Sec. IXII we developed a 
new approach to the description of neutrino oscillations 
in vacuum which is based on the consideration of fla- 
vor neutrinos emitted by classical sources. Within this 
method one could study neutrino oscillations in space 
of both Dirac [Eqs. (I11.24j) and (|11.25|) ] and Majorana 
[Eqs. (|11.46[) ] neutrinos. For the case of Majorana neu- 
trinos we could also describe the transitions between neu- 
trinos and antineutrinos [see Eq. (| 1 1 .4T[) ] . This our result 
is also consistent with the previous studies [2g, |92[ . In 
frames of this formalism we exactly take into account 
neutrino masses as well as examine what kind of neu- 
trino sources could give us a stable oscillations picture 
(see also Appendix ICj) . 

However, the external classical source approach is not 
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free from disadvantages. This method is based on the use 
of the retarded Green functions [Eqs. (|11.8|) and (|11.34[) ] 
in the neutrino mass eigenstates basis. Thus, if one tries 
to apply this formalism for the description of neutrino 
oscillations in an external field, which, as we mentioned 
above, is typically nondiagonal in the mass eigenstate 
basis, one encounters a problem to find a nondiagonal 
retarded Green function. It is quite difficult to calculate 
this kind of function analytically. Thus this approach 
is unlikely to be applicable to the description of neu- 
trino oscillations in an external field, in contrast to the 
relativistic quantum mechanics method with gives satis- 
factory results in presence of external fields at least for 
ultrarelativistic neutrinos, unless we study a special case 
of an external field which is diagonal in the neutrino mass 
eigenstates basis. 

We made an attempt to describe transitions between 
neutrinos and antineutrinos using the relativistic quan- 
tum mechanics method [Eq. (|8.13l) ] and the external clas- 
sical source approach [Eq. (jll.47j) ] and obtained the cor- 
responding transition probabilityformulae which resem- 
ble the previously derived ones [26|,|92|. However, in case 
of Majorana neutrinos, particles are identical to their an- 
tiparticles. The concept of antiparticles is an inherent 
component of quantum field theory. Thus the approaches 
based on classical physics do not seem to be an appropri- 
ate tool for the description of oscillations between neu- 
trinos and antineutrinos. That is why in Sec. IXHI we 
use the quantum field theory method to study particles 
and antiparticles in background matter in case of Majo- 
rana neutrinos. Note that the analogous method for the 
description of neutrino oscillations in vacuum was previ- 
ously used in Refs. [HHzlIil. 

It should be noted that oscillations between neutrinos 
and antineutinos manifest in (Oi'2/3)-decay 17]. During 



this process two nucleons inside a nucleus exchange a 
virtual Majorana neutrino. That is why it is important 
to examine the influence of dense nuclear matter on the 
process of Majorana neutrinos propagation. 

In our analysis we considered the situation when the in- 
teraction with matter is diagonal in the mass eigenstates 
basis since the case of a nondiagonal neutrino propaga- 
tor is rather difficult to study analytically. This diffi- 
culty is analogous to that in the external classical source 
approach. Nevertheless we found that this type of mat- 
ter interaction suppresses transitions between neutrinos 
and antineutrinos. Moreover, in case of hypothetical very 
heavy Majorana neutrinos the transition probability van- 
ishes [see Eqs. (|12.11[) and (|12.12[) ]. It means that one 
cannot explore the presence of this kind of heavy neutri- 
nos studying (0^2/3)-decay. 

In Sees. IVI1 IVII1 and [X] we discussed several applica- 
tions of the results obtained in frames of the relativistic 
quantum mechanics method to the studies of astrophys- 
ical neutrinos emitted during the core collapse of a su- 
pernova. First we studied spin flavor oscillations of Dirac 
fSec. lVI[) and Majorana (Sec.|X|) neutrinos in an expand- 
ing envelope under the influence of the magnetic field of 



a supernova. It was found that for Dirac neutrinos the 
amplification of neutrino oscillations can be achieved in 
a moderate magnetic field [see Figs. [Ud)-|3Jd)]. On the 
contrary, in the Majorana neutrinos case, neutrino oscil- 
lations can be enhanced only in a strong magnetic field 
[see Fig. OP), (d) and (f)]. 

In Sec. Eni another channel of oscillations of astrophys- 
ical neutrinos was considered. We studied the possibil- 
ity of spin flavor oscillations between electron and ad- 
ditional quasi-degenerate in mass sterile neutrinos. We 
obtained that the correction f|5 . 1 7[) to the effective Hamil- 
tonian (I5.16[) , found in frames of the relativistic quantum 
mechanics method, results in the appearance of the new 
resonance in this oscillations channel. 

Finally in Appendix [X] we present the general solution 
of the ordinary differential equations for the coefficients 
which one encounters in the relativistic quantum me- 
chanics method. In Appendix [B] we discuss the validity 
of the relativistic quantum mechanics approach to the 
description of neutrino spin flavor oscillations. In par- 
ticular, we analyze various factors which influence oscil- 
lations between electron and sterile neutrinos with very 
small 8m 2 (see also Sec. IVII| ). In Appendix [Cl we study 
how neutrino wave functions converge to the values cor- 
responding to the stable oscillations picture in the exter- 
nal classical source approach (Sec. IXI|) . In Appendix [D] 
we present the details of the calculation of the S'-matrix 
element involving transitions between neutrinos and an- 
tineutrinos in frames of the quantum field theory descrip- 
tion of neutrino oscillations (Sec. IXlIj) . 

Summarizing we can say that a theoretical approach 
for the description of neutrino oscillations should satisfy 
the following requirements: the evolution equation gov- 
erning the dynamics of mixed neutrinos should be de- 
rived from Lorentz invariant Lagrangian of the system 
and thus it should account for the relativistic invariancc, 
at least implicitly; and such a method should be applica- 
ble in presence of various external fields. Among various 
theoretical approaches [25|, [96| the relativistic quantum 
mechanics method is likely to be most appropriate for- 
malism for the description of neutrino oscillations. 
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Appendix A: Solution to the ordinary differential 
equations for the functions ai^ 

In this Appendix we describe the formalism for the 
analysis of ordinary differential equations for the func- 
tions a„ which one encounters in Sees. Illll and IIV1 
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Let us study the time evolution of the two compo- 
nent spinor Z T = (Zi,Z 2 ) which is governed by the 
Schrodinger equation of the form, 

iZ = HZ, (Al) 

where the Hamiltonian has the following form: 



H = g 











(A2) 



Here g and ui are real parameters. Eq. (|A1[) should be 
supplied with the initial condition Z(0). To find the so- 
lution of Eqs. (|A1[) and (|A2I) we introduce the new spinor 
Z' by the relation, Z = UZ' , where the unitary matrix U 
reads 



U 



' e iwt/2 q 

e" iwt / 2 



(A3) 



Now Eq. (|Alj) is rewritten in the following way: 

iZ' = H'Z', (A4) 

with the new Hamiltonian H ' which is obtained with help 
of Eqs. JXTJ-dAD, 



H' = U^HU- SUtU = 




(A5) 



Note that the initial condition for the spinor Z'(0) is the 
same as for Z(0), Z'(0) = Z(0), due to the special form 
of the matrix U in Eq. ([A3]) . 

Supposing that the Hamiltonian H' in Eqs. (IA4[) 
and (IA5|) does not depend on time we get the solution to 
Eq. ([All) as 

Z'(t) =exp(-Lff't)Z'(0) 

= (cos m - i(o-n) sin J?f) Z'(0), (A6) 

where n = (g,0,u>/2)/ ft is the unit vector and ft = 
yjg 2 + (w/2) 2 . Using Eqs. (jA3]l and (|A6]l we arrive to 
the expressions for the components of Z written in terms 
of the initial condition Z(0): 

Z 1 (t) = (cos fit - i^j sin J?t) e iwt/2 Zi (0) 



i— sm^e^* 72 ^), 



. 9 ■ 

Z 2 (t) = (cos Qt + sin fit) e _i "' /2 Z 2 (0) 
-i^sin(m)e- iut / 2 Z 1 (0), 



(A7) 



which can be directly applied for the analysis of ordinary 
differential equations from Sees. Illll and IIVI 

To get the solution of Eq. (|3.11[> we identify the compo- 
nents of the spinor Z with a^ 2 and the parameter u> with 
<jj- (see Sec. IIIip . Finally we arrive to Eq. (|3.12|) . We can 
also apply Eq. (|A7I) to obtain the solution of Eq. (|4.12l) . 
For this purpose one considers two cases: 



• For Z T 



(af,a£) 



and fi 



• For Z T 



'1 i "2 



a 2 ), g = fJ-B, uj = uj- and fl = f2_. 



Using these formulae together with Eqs. (|A7[) we readily 
arrive to Eqs. (|4.13j) - (|4.15j) . Note that the dynamics of 
the system (|Alj) and (|A2|) is analogous to the quantum 
mechanical description of neutrino spin flavor oscillation 
in a twisting magnetic field studied in Ref. |97| . 



Appendix B: Analysis of approximations made in 
the derivation of the effective Hamiltonian 



In this Appendix we analyze the validity of the rela- 
tivistic quantum mechanics approach for the description 
of neutrino spin flavor oscillations used in Sec. [Vj The 
correction to the effective Hamiltonian (|5.17|) is a rather 
small quantity. Therefore we should evaluate other fac- 
tors which can also give the contributions, comparable 
with Eq. (|5.17p , to the effective Hamiltonian. In this sec- 
tion we analyze the contributions to the effective Hamil- 
tonian from longitudinal magnetic field, matter polariza- 
tion and possible corrections from the new interactions. 

First we should remind that we use the relativistic 
quantum mechanics approach, with the external fields be- 
ing independent of spatial coordinates. If external fields 
depend on the spatial coordinates, Dirac wave packets 
theory reveals various additional phenomena such as par- 
ticles creation by the external field inhomogeneity 98| . 
For the approximation of the spatially constant exter- 
nal fields to be valid, the typical length scale of the ex- 
ternal field variation L cx t should be much greater than 
the Compton length of a neutrino [98j |: L cx t %c = 
H/m^c [49(. For a neutrino with m v ~ leV this condi- 
tion reads L ext 10~ 5 cm, that is fulfilled for almost all 
realistic external fields. 

A general remark on the "perturbative" approach used 
in Sec. IVl should be made. In the modified Eq. Q3.9p . 
which includes the interaction with the magnetic field, 
the terms containing and aj are coupled and the cou- 
pling terms are proportional to g and fiB (|5.8p . As we 
demonstrate in Sec. HH mass eigenstates decouple in vac- 
uum are their values depend on the initial condition only. 
While solving the modified Eq. (|3.9p . we could just take 
a term linear in g and /j,B (see, e.g., Eq. (|4.20l) as well as 
Refs. [13, HI])- However we take into account these terms 
exactly in the further analysis [see Eqs. (|5.7p and (|5.10p ]. 
It is equivalent to the summation of all terms in the per- 
turbation series. 

While deriving the effective Hamiltonian (|5.I0I) in 
Sec. [V] we supposed that magnetic field is transverse with 
respect to the neutrino motion. The effect of the longitu- 
dinal magnetic field on neutrino oscillations was studied 
in Ref. [55|. It was found there that diagonal entries of 
the effective Hamiltonian receive additional small contri- 
butions [i a B\\ (m a /k). In order to neglect the longitudinal 
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magnetic field contribution in comparison with our cor- 
rections (I5.17[) , its strength should satisfy the condition, 



Bn 



1 



B±_ l6kB±\[i a m a - [i b m b \ 



„,2 3 
m a9a 



m b9b 



Ml 



Ml 



(Bl) 



where B± is the transverse component of the magnetic 
field. 

A regular magnetic field of a supernova typically has 
both poloidal and toriodal components. Of course, an 
irregular turbulent magnetic field can be also present but 
its length scale appears to be quite small. A toroidal 
magnetic field can be ~ 10 16 G and is concentrated near 
the equator of the star at the distance ~ 1 km from the 
star center [99] . Thus in our case a toroidal magnetic field 
is unlikely to significantly contribute to the dynamics of 
neutrino oscillations. 

Let us evaluate the fraction of neutrinos for which the 
new correction (|5.17[) to the effective Hamiltonian gives 
bigger contribution to the resonance enhancement of os- 
cillations compared to that of the longitudinal compo- 
nent of the poloidal magnetic field. Using Eq. (|B1[) we 
find that these neutrinos should be emitted inside the 
solid angle near the equatorial plane with the spread 2§, 
where i? ~ Bn/B_\_- Assuming the radially symmetric 
neutrino emission we find that about 10% of the total 
neutrino flux is affected by the new resonance (I7.1[) . i.e. 
the influence of the longitudinal magnetic field is negli- 
gible for oscillations of such particles. 

The next important approximation made in the de- 
viation of Eq. (15 . 10[) was the assumption of negligible 
polarization of matter which can be not true if we study 
rather strong magnetic fields. The effect of matter polar- 
ization on neutr ino o scillations was previously discussed 
in Refs. (49l. [Hoi . flOOj ] . Matter polarization produces the 
following contributions to the diagonal entries of the ef- 
fective Hamiltonian [49|, <7a(A//3„) (left polarized 
neutrinos) and <?a(A//3„)(m /fc) (right polarized neutri- 
nos) . Here (3 V is the neutrino velocity and A f is the mean 
polarization vector of background fcrmions and we keep 
only the leading order in m a /k. 

First we estimate the contribution to the right polar- 
ized neutrinos effective potential. It is clear that one 
should take into account only the polarization of elec- 
trons since nucleons are much heavier. For the weakly 
degenerate electrons we s hould discuss the case of weak 



10" 



< 1, 



field limit (see Ref. (lOOj), since 2eB/m 
where m e is the electron's mass and B ~ 10 7 G. Since the 
temperature i nsid e the shock wave region can be about 
several MeV jl0l| the electrons are relativistic. Hence 
their mean polarization can be estimated as, |A/| ~ 
/iB-Bm e /3Tg , where T e is the temperature of electrons. 
Therefore we get that the new correction to the effective 
Hamiltonian (|5 . 1 7[) becomes bigger than the contribution 
of matter polarization to the effective potential of right 
polarized neutrinos at T e > 4 MeV. 

To evaluate the contribution g a (\ff3 v ) to the effective 
Hamiltonian, which corresponds to the left polarized neu- 
trinos effective potential, we notice that the vector A/ 



should be directed along the external magnetic field. For 
this term to be much less than the new correction to the 
effective Hamiltonian (|5.17p . the angle defined above, 
should be very small: d < 10 -8 , for T e ~ 4 MeV. It 
means that we can neglect the polarization effects only 
if neutrinos are emitted very close the equator of a star. 
We should, however, remind that this kind of matter po- 
larization term contributes only to {Hqm)22 in Eq. (|5.16[) 
since only this entry corresponds to v\. Thus the pres- 
ence of the term g a (\f(3 u ) does not directly affect our 
results since we study vf- o oscillations channel (see 
also Table H|. 

The presence of big Dirac neutrino magnetic moments 
implies the existence of new interactions, beyond the 
standard model, which electromagnetically couple left 
and right polarized neutrinos. It is probable that these 
new interactions also contribute to the effective poten- 
tial of the right polarized neutrino interaction with back- 
ground matter. Despite this additional effective potential 
is likely to be small, one should evaluate it and compare 
with the correction (|5.17|) . 

The most generic SU(2)l x U(l)y gauge invariant and 
remormalizable interaction which produces Dirac neu- 
trino magnetic moment was discussed in Ref. [l02j |. It 
was found that neutrino magnetic moments arise from 
the effective Lagrangian involving the dimension n = 6 
operators Oj, 



h.c. 



(B2) 



where A ~ 1 TeV is a scale of the new physics and Cj are 
the effective operator coupling constants. The sum in 
Eq. (|B2[) spans all the operators of the given dimension. 

One of the operators Oj also contributes to the effec- 
tive potential of a right-handed neutrino in matter, 



(B3) 



d^W? - d u W« 



where k is the coupling constant, r a are Pauli matri- 
ces, L T = (z/jj , &l) is the S\J(2) L isodoublet, <f> = iT 2 </>*, 
with <fi being a Higgs field, and W^ v 
ne abc W*Wc is the SU(2) L field stren 

Assuming the spontaneous symmetry breaking at the 
electroweak scale, 4> T —> 
Lagrangian in the form, 



(0,v/^/2), we can rewrite the 



Ckv 



off 



(B4) 



v\j should 



which implies that a process e~ + j^r — > e" 
happen in background mat ter. 

Using the results of Ref. [l02j we can evaluate the con- 
tribution of the Lagrangian (|B4[) to the effective Hamil- 
tonian (|5.16[) as 



6V R ~ V„ 



Mb 



(B5) 
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discuss two functions 




x, 10° km 



FIG. 7. The numerical transition probability for neutrinos 
with p, = 3x 10 _1 Vb, /it„ e = 10~ 13 mb, E = 100 MeV, m Ve = 
leV and Sm 2 — 5 x 10 _17 eV 2 moving in background matter 
with p = 10 6 g/cm 3 and interacting with the magnetic field 
B = 10 7 G. After Ref. \M- 



where V sm ~ Gpn e is the standard model effective poten- 
tial and Mw is the W boson mass. Taking p v ~ 10~ 12 ^b, 
E v ~ 100 MeV and m v ~ 1 eV (see Sec. |VTI|) we can get 
that the ratio of the correction to the effective poten- 
tial (|B5j) and new correction (|5.17|) is ~ 10~ 4 . It means 
that the influence of new interactions, which generate 
neutrino magnetic moments, are not important for neu- 
trino spin flavor oscillations. 

Now let us estimate the influence of the diagonal mag- 
netic moments p a on the dynamics of spin flavor oscil- 
lations. It was found in Ref. [(36[ that to get the signif- 
icant luminosity ~ 10 50 erg/s the diagonal magnetic 
moment should be fx Ve = 10~ 13 ^b- Eq. (|7.1j) was ob- 
tained under the assumption /i„ e -C p. In Eq. (|7.3p we 
use /i = 3x 10 _12 ^b, i-e. the condition of the Eq. (|7.1[) 
validity is satisfied. In Fig. [7] we present the numerical 
solution of the Schrodinger equation (|5 . 10[) . We remind 
that our simplified model with /i„ e = gives the transi- 
tion probability P(x) — sin 2 (pBx) if the resonance con- 
dition (17.11) is fulfilled. As one can see on Fig. [7] there is 
almost no difference in the dynamics of spin flavor oscilla- 
tions in our simplified model and more realistic situation 
which involves non-zero diagonal magnetic moment of an 
electron neutrino. 



Appendix C: Evaluation of integrals 

In this Appendix we calculate the integrals which one 
encounters while studying neutrino oscillations in the 
model with classical sources in Sec. IXII It is interest- 
ing to evaluate the inexactitude which is made when we 
approach to the limit x m — ¥ oo in Eq. (|11.15[) . Let us 



Vr 1 



"2/ 



F Q (r) =-(e i£ 
P 



(CI) 
(C2) 



These functions are proportional to I2 in Eqs. (|11.15|) 
and pi,16p respectively. In Eqs. (|C1|) and (IC2[) we use 
dimensionless parameters p = mr, £ = Er = yp, 7 = 
E/m, y m = y/(t/r)* - 1 and P = V^ 2 - P 2 ■ 

In case we study neutrinos, we get that £ ^> p ^> 1 
in almost all realistic situations. For example, suppose 
we study a neutrino emitted in a supernova explosion is 
our Galaxy. The typical distance is r ~ lOkpc. Taking 
m ~ 1 eV and E ~ 10 MeV, we receive that £ ~ 10 35 
and p ~ 10 28 . 

Basing on the analysis of Sec. IXI Al we can rewrite 
Eq. JUB) as 



F(r,t) =F (r)-6F, 
5F = I dy 



(C3) 



Using the fact that p> 1, y m 3> 1 and the representation 
for the Bessel function, 



3tt 

T 



at 



-00, 



(C4) 



we obtain for the function SF in Eq. (|C3[) the following 
expression: 

SF -^={ci([j + l]py m ) + ci([7 - l]py m ) 



+ i[si([ 7 + l]py m ) + si([ 7 - l]py m )}}, (C5) 

where ci(z) and si(z) are cosine and sine integrals. Using 
the asymptotic expression, 

|ci(z)| ~ |si(z)| ~ - at z -> +00, (C6) 

z 

we obtain that the function 5F approaches to zero as 
l/(y m p 3 / 2 ) at great values of y m and p. Note that this 
result remains valid for a particle with an arbitrary 7 fac- 
tor, i.e. rapid oscillations of the function SF will atten- 
uate even for slow particles. This analysis substantiates 
the approximations made in Sec. IXI Al 

Finally let us illustrate the behavior of the functions 
F(r,t) and Fo(r). On Fig. [8] we present the absolute 
values of these functions versus t. This figure is plotted 
for p = 1 and £ = 10. The solid line is the absolute 
value of the function F(r,t) and the dashed line of the 
function Fo(r). As we mention in Sec. IXI Al thc relaxation 
phenomena occur when t > r. It can be also seen on 
Fig. [5] It is possible to notice that \F(r, t)\ — > \F (r)\ at 
great values of t as it is predicted in Sec. IXI Al 
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FIG. 8. The absolute values of the functions F(r, t) and Fo(r) 
versus t. This figure is taken from Ref. [3||. 



Appendix D: Calculation of the S- matrix element 

In this Appendix we present the detailed calculation of 
the .S-matrix element in Sec. IXIII Note that analogous 



calculation for the case of virtual neutrinos propagat- 
ing in vacuum, i.e. when the matrix is absent in 
Eq. (fT2~T]) . is presented in Ref. [2rjj ]. 

Using Eqs. (|12.5[) - (|12.T[) as well as after the spatial in- 
tegration and the elimination of the combinatorial factor 
we can cast Eq. (|12.3|) in the form, 



S=- 



dx dy e 



iP/3Xi-ip Q x 2 iE a x -iEfjy 



x u( Pa h P R S a (L,x - y )PR7 Q u( Pp ), (Dl) 



where we use the fact that u T (— p a )C — u(p a ). The 
integration with respect to xq and yo can be performed 
with help of the new variables: T = (xq + yo)/2 and 
t = x — y . After this integration one gets the energy 
conservation ^-function in Eq. (|12.8[) . 

The Fourier transform of the ne utrin o propagator 
S a (x) in Eq. (|Dlj) was found in Ref. |l03| and has the 
form, 



S a (k) = 



(fc 2 - ml - g|/4 - ky^ff^)^ + m + 7„7 5 ff£/2) 
(fc 2 - m 2 - <? 2 /4) 2 - (<? a fc) 2 + k 2 g 2 a 



(D2) 



r 



where g% = g£ a are the diagonal elements of the matrix 
04) [seeEq. ([3D]. 

In case of the nonmoving and unpolarized matter the 
momentum integration in Eq. (|D1[) can be performed us- 
ing the calculus of residues. Finally we arrive to the 
following result: 



r .kL d 3 k pi - gl/4 - k 2 - g a (ak) _ . 



(2tt) 3 (k 2 — fc 2 )(k 2 - fe| 



(D3) 



where i 7 ^ is defined in Eq. (112. 9p . In Eq. ()D3|) we neglect 
several small terms since k\^L 1. 

It is convenient to perform momentum integration in 
Eq. (|D3[) using cylindrical coordinates pointing L along 
the z-axis. Note that the integration over fc^ is trivial 
and gives 2-7T since the terms in the integrand [see also 
Eq. (|D2[) ] containing fc^ just vanish. The result of the 



integration over k can be presented in the form, 
d 2 k F(k pi k z ) 

(2tt) 2 (k 2 - fc 2 )(k 2 - fc 2 ) ~ a + b+ c ' [ ' 
where d 2 k = k p dk p dk z and 

j — I 2 kpdkp l' + ^ F(kp, k z ) 



(27T) 2 



: lh.2 _ U2 \(U2 _ U2 \ ' 
\^z ^zljK^z K z2) 



fel 



fe2 (27T) 2 7^ 



■/, = / fe^fe ^dfc ^ F L^ z) (D5) 



(fc 2 + fc" 2 )(fc 2 + fc"|) 



fel (2^) 2 



Here we use the notations fc^ 2 = 1/^12"^ an d 2 = 



fc 2 — fc 2 2 . The contours of the integration for each of 
the integrals in Eq. (|D5I) are shown in Fig [SJ 
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